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ABSTRACT 

We present the results of Spectral Energy Distribution (SED) fitting analysis for Lyman Break 
Galaxies (LBGs) at z - 5 in the GOODS-N and its flanking fields (the GOODS-FF). With the 
pubhcly available IRAC images in the GOODS-N and IRAC data in the GOODS-FF, we constructed 
the rest-frame UV to optical SEDs for a large sample (~ 100) of UV-selected galaxies at z 5. 
Comparing the observed SEDs with model SEDs generated with a population synthesis code, we 
derived a best-fit set of parameters (stellar mass, age, color excess, and star formation rate) for each 
of sample LBGs. The derived stellar masses range from 10^ to IO^^Mq with a median value of 
4.1 X 1O^M0. The comparison with z = 2 — 3 LBGs shows that the stellar masses of z ~ 5 LBGs 
are systematically smaller by a factor of 3 — 4 than those of z — 2 — 3 LBGs in a similar rest-frame 
UV luminosity range. The star formation ages are relatively younger than those of the z = 2 — 3 
LBGs. We also compared the results for our sample with other studies for the z — 5 — 6 galaxies. 
Although there seem to be similarities and differences in the properties, we could not conclude its 
significance. We also derived a stellar mass function of our sample by correcting for incompletenesses. 
Although the number densities in the massive end are comparable to the theoretical predictions from 
semi-analytic models involving AGN feedback, the number densities in the low-mass part are smaller 
than the model predictions. By integrating the stellar mass function down to 10^ Mq, the stellar mass 
density at z ~ 5 is calculated to be (0.7 — 2.4) x lO^MoMpc"^. The stellar mass density at z ~ 5 is 
dominated by massive part of the stellar mass function. Compared with other observational studies 
and the model predictions, the mass density of our sample is consistent with general trend of the 
increase of the stellar mass density with time. 

Subject headings: galaxies: high-redshift — galaxies: evolution — galaxies: stellar content — galaxies: 
luminosity function, mass function 



1. INTRODUCTION 

A large number of Lyman Break Galaxies (LBGs) at 
z = 2 — 4 have been found and studies of their prop- 
erties have been made extensively fe.g .. ISteidel et al] 
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time have been revealed. Particularly, the evolution 
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of the UV-luminosity function and the UV-luminosity 
density of galaxies from z 6 to present are extensively 
studied. Since the rest-frame UV luminosity density 
traces star-formation activity in galaxies, the recent 
studies are revealing the cosmic sta r-formation history. 
A com pilation of the results by iHopkins fc BeacomI 
()2006f ) shows that the cosmic star formation rate (SFR) 
density increases from z = 6 — 10 to z = 2 — 3 and 
decreases to present. It is also important to investigate 
the stellar mass assembly in the high-redshift universe. 
It is thought that galaxies are assembled from smaller 
systems through mergers in the framework of the cold 
dark matter (CDM) hierarchical structure formation 
scenario. Thus a stellar mass of an individual galaxy 
and a stellar mass density of galaxies are expected to be 
increasing with redshift. The stellar mass together with 
its age, color excess, and star formation rate are esti- 
mated by fitting of spectral energy distribution (SED) 
of an LBG (Sawicki fc Yee 19981. This method is now 
widely used t o constrain stel l ar p opulati on of galaxies. 
For in stance, iPapovich et al.l (jlpOl) and IShaplev et al.l 
(|2001f ) studied stellar population of LBGs at z ~ 3 and 
found that the LBGs have stellar masses of ~ IO^^Mq 
and they are dominated by relatively young (several 
tens to several hundreds of Myr) stellar populations. 
Similar study wa s made at z ~ 2 (^ 1 Gyr after z ~ 3) 
for BX galaxies (jShaplev et al.l [20051 ): the massive side 
of stellar mass distribution at z 2 seems to increase 
slightly, that might indicate the stellar mass evolution 
of an LBG. Stellar masses for other kinds of galaxy 
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populations at z = 2 — 4 as well as for galaxy samples 
with photometric redshifts have also been studied; these 
observations show the gradual growth of th e stellar mass 
density at z < 5 (e .g. Shapley et al.l 120051 iDaddi et al. 
20041: iFra nx et al.' '2003'; 'Brammer & van Dokkum 



20071 iBrinchmann fc EUis.2000; Dickinson et al. 2003b: 



Drorv et all 120051: iFontana et al.l l2006t iRudnick et all 



2006HPerez-Gonzalez et al.ll2008t fsisner et al.ll2008[) . 

Most of the stellar masses of z < 5 galaxies have been 
derived based on observations at wavelengths from op- 
tical to near- infrared (NIR); NIR data is necessary to 
cover the rest-frame wavelength region larger than 4000 
A, which is sensitive to the stellar mass of a galaxy. 
In order to push the redshift higher, the observations 
in mid-infrared are required to constrain stellar masses. 
With the advent of the Spitzer Space Telescope (hereafter 
the Spitzer) we can now acc ess the longer wa velengths. 
The Infrared Array Camera (|Fazio et al.ll2004 hereafter 
IRAC) on the Spitzer allows us to reach rest-frame op- 
tical wavelength for z > 5 galaxies. Only recently, there 
have been studies of the stellar populations of galaxies 
at z > 5 in the southern field of the Great Obse rvato- 
ries Origins Deep Survey ([Dickinson et al.ll2003 al. here- 
after the GOODS-S), where deep optical to MIR images 
obtained w ith IR AC are av ailable including NIR data 
(lYan et all 2006; Labbe et al. 2006; Stark ct al. 2003; 
iVerma et aO 120071: lEvles et al.l 12007). However, the re- 
sults are only obtained in the GOODS-S. Considering 
relatively small sizes of the samples and a possible cos- 
mic variance, we need to study with a larger sample and 
in another field to reach robust consensus. 

In this paper, we reveal stellar populations, espe- 
cially stellar mass, of LBGs at z '--^ 5 in a field con- 
taining the Great Observator ies Origins Deep Survey 
North (jDickinson et all l2003al . hereafter the GOODS- 
N), where our group constructed a large sample of LBGs 
at z ~ 5, studied the eyoluti on of UV- luminosity func- 
tion (jlwata et al.ll2003l I2007D. and made spec troscopic 
follow-up observations (|Ando et al.ll2004 [20071 ). In ad- 
dition to public IRAC data in the GOODS-N, we made 
Spitzer/IRAC observations in the flanking fields (here- 
after the GOODS-FF) of the GOODS-N. Although the 
depth of the GOODS-FF is 1 ~ 1.5 mag shallower than 
the GOODS-N, it gives us more than twice as large area 
as the GOODS-S and consequently making the number 
of our sample galaxies large. This would especially be 
important to increase the sample size of massive side of 
galaxy stellar mass function. We perform the SED fit- 
ting analysis with this large sample and derive the stellar 
populations of galaxies and aim at investigating the evo- 
lution of stellar populations from z~5toz = 2 — 3. 
We also derive the stellar mass function and density at 
z ~ 5. This paper is organized as follows. In §2 we 
summarize data and sample selection. We outline our 
population synthesis modeling in §3 and present results 
in §4. In §5, comparison with LBGs at lower and higher 
redshifts and estimation of stellar mass function and stel- 
lar mass density are presented. Summary is given in §6. 
Throughout this paper, we adopt the concordance cos- 
mology, (Ojvf, ^ A, h) = (0.3, 0.7, 0.7 ). All magnitudes are 
on AB system ([Oke fc Gunn|ll983f) . 

2. DATA AND LYMAN BREAK GALAXY SAMPLE 



2.1. Optical Data 

The opti cal imaging data used in this work are 
taken from llwata et al.l ([2007f) . They carried out deep 
and wide imaging observations of two independent 
blank fields, namely, the field including the GOODS- 
N and the J0053-H 234 field, with the Suprime-Cam 
([Mivazaki et al.|[2002f) attached to the Subaru Telescope. 
Here we use a sub-sample of the GOODS-N field with 
an effective survey area of 508.5 arcmin^. The im- 
ages were taken through the V, Ic, and z' filters. The 
FWHMs of the reduced data are ~ 1.1". The detailed 
descriptions of the ima ge properties and t he data reduc- 
tion are presented bv llwata et al.l ([20071 ). Object de- 
tection and photom etry were made by using SExtractor 
([Bertin fc Arnoutd [r996). They measured MAG_AUTO 
and 1.6" diameter aperture magnitudes for total mag- 
nitudes in z'-band and Ic — z' colors, respectively. For 
total magnitudes in /c-band, we calculated them from 



total magnitudes in z'-band and Ic 



colors. V^-band 



magnitudes are obtained by 1.6" diameter aperture pho- 
tometry. The limiting magnitudes for V, 1^ and z'-band 
images are 28.1, 26.8, and 26.6 mag, respectively (3ct, 
1.6" diameter aperture). We do not use the T^-band mag- 
nitudes in the SED fitting analysis. 

2.2. Mid- Infrared Data 

A part of the Subaru imaging area was observed with 
the IRAC. We use the pubhcly available GOODS-N 
IRAC data in this work. We use the First Data Release 
(DRl) and the Second Data Release (DR2) of IRAC data 
products from the GOODS Spitzer Legacy Science pro- 
gram, which consist of imaging data in 3.6, 4.5, 5.8, and 
8.0 fim passbands with the total effective area of ~150 
arcmin^. The 3a limiting magnitudes in 2.4" diameter 
apertures are 25.9, 25.6, 23.7, and 23.6 mag in the 3.6, 
4.5, 5.8, and 8.0 ^m bands, respectively. 

In addition to the GOODS-N data, we obtained IRAC 
data for the GOODS-FF to cover the most part of 
the Subaru imaging area. The IRAC data in the the 
GOODS-FF were obtained in 2005 December and 2006 
June as General Observer (GO) program 20218, and are 
1 — 1.5 mag shallower but ~100 arcmin^ wider than the 
GOODS-N IRAC imaging data. 

We used the Basic Calibrated Data processed by the 
IRAC data reduction pipelines (version S14.0.0) at the 
Spitzer Science Center (SSC), and the MOPEX pack- 
age (version 030106) was used for further reduction. Af- 
ter removing artifacts (mux bleed and column pulldown) 
on the images, background counts of individual images 
are estimated and subtracted. Then the pointing refine- 
ment was made to improve the consistency of positions 
of individual images. Finally, the individual frames were 
drizzled and mosaiced to create a single image of the 
GOODS-FF for each IRAC channel. The pixel scale of 
the mosaiced images was set to be 0.606", which is ap- 
proximately half of the native IRAC pixel scale and al- 
most comparable to the pixel scale of the public images 
in the GOODS-N. The 3cr hmiting magnitudes in 2.4" 
diameter apertures are 24.8, 24.1, 22.2, and 22.3 mag in 
the 3.6, 4.5, 5.8, and 8.0 fim bands, respectively. Com- 
bining the GOODS-N IRAC data and the GOODS-FF 
data, we covered a total effective area of ^400 arcmin^, 
which covers ~80% of the area taken with Subaru. 
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Fig. 1. — Postage stamps (20"x20") of three representative 
objects in 5 passbands. From left to right, V, Ic, z' , 3.6/xm, and 
4.5/im-band images are shown and the LBG candidate is indicated 
by a cross in each panel. North is up and east is to the left. 

Source detection and photometry were made by us- 
ing SExtractor. As we discuss below, because the IRAC 
images are crowded and neighboring objects may affect 
photometry, we performed aperture photometry in aU 
channels with a diameter of 2.4" and applied aperture 
corrections to obtain total magnitudes. We examined 
the best value of the aperture size and chose the value 
of 2.4", which maximizes the signal-to- noise ratio(S/N). 
The correction factors from aperture magnitudes to total 
magnitudes were derived from Monte Carlo simulations 
in which artificial objects with IRAC point spread func- 
tions (PSFs) were put into the images, then detected 
and measured with the same SExtractor parameters as 
we adopted. In these simulations, our targets were as- 
sumed to be point sources because their apparent size is 
small enough as compared to the PSF. The PSF of the 
IRAC images was made by stacking the IRAC images 
of objects showing SExtractor "Stellarity" indices (1 for 
point sources and for extended sources) larger than 0.98 
in the Subaru data. For the sample of the GOODS-N, 
the factors are —0.69, —0.72, —0.99, and —1.06 mag in 
3.6, 4.5, 5.8, and 8.0 /im bands, respectively. The uncer- 
tainties of the correction are ~ 2%, ~ 3%, ~ 10%, and 
~ 10% in 3.6, 4.5, 5.8, and 8.0 /zm bands, respectively. 
For the sample of the GOODS-FF, the factors are -0.70, 
-0.73, -0.95, and -1.01 mag in 3.6, 4.5, 5.8, and 8.0 fiui 
bands, respectively. The uncertainties of the correction 
are ~ 3%, ~ 5%, ~ 30%, and ~ 30% in 3.6, 4.5, 5.8, and 
8.0 /Ltm bands, respectively. The errors of the resulting 
magnitudes are taken to be la standard deviation of sky 
background. Because of the too low S/N to provide use- 
ful upper limits in the SED fitting and large uncertainty 
of the correction factors in the 5.8/im and 8.0/im bands, 
we do not use data of these bands for the SED fitting. 

2.3. Sample 

The LBG sample we use in this work is selected b y the 
following color selection criteria (jlwata et al.ll2007[ ): 

V -Ic> 1.55 & V -Ic> 7.0(4 - z') + 0.15 . (1) 

The number of objects which satisfy these selection cri- 
teria is 617 in z' < 26.5 mag. These criteria are de- 
signed to select LBG candidates at z ~ 5 (hereafter, we 
refer them as LBGs at z ~ 5, though they are candi- 
dates in the strict sense.) efficiently without heavy con- 
tamination by interlopers such as objects at lower red- 
shift and stars in the Galaxy. Follow-up spectroscopy of 



the candidates confirms that the selection c riteria effec- 
tively extra ct star-forming galaxies at z ~ 5 (jAndo et al.l 
|2004 120071 ). though, the number of spect roscopically 
confirmed o bjects is still limited to ^ 10. Ilwata et al.l 
(|2003l l2007| ) estimated a fraction of the interlopers em- 
ploying a resampling method, and found the estimated 
fractions of interlopers are ~ 50%, ~ 20%, and ^ 10% in 
z' = 23.0 - 24.0 mag, 24.0 - 26. mag, and 26.0 - 26.5 
mag, respectively (see Table 5. of Ilwata et al.ll200"7h . We 
correct for the factors, when we derive the stellar mass 
function in 

There are dee p X-ray observations with Chandra in 
the GOODS-N (|Alexander et al.l[200l . About 60% of 
our sample LBGs lie in the region covered by Chan- 
dra. We cross-matched the LBGs with X-ray sources and 
found one object is an active galactic nucleus (AGN) at 
z = 5.186 with Lx = 6.8 x 10^^ g-i (2-10 keV) 
(|Barger et al.l 120021 : lAndo et al.ll2004[ ). We do not use 
this object in the SED fitting analysis. All others are not 
detected at the 3cr flux limit of 2.5 x 10~^^ erg s""'^ cm~^ 
(0.5-2.0 keV), which corresponds to Lx 6 x 10^^ erg 
s~^ (2-10 keV). This luminosity level is that for Seyfert 
class AGNs, and hence the LBGs do not harbo r X-ray lu- 
minou s AGNs or they may be obscured AGNs. iHasingerl 
(2008') suggested that the fraction of type-2 AGN in- 
creases with redshift, and at z = 3 — 5 the type-2 fraction 
in luminosity range of Lx = 10^^ — lO''^ erg s~^ is ~ 0.9. 
Hence, in any cases, AGN components presumably do 
not affect the SEDs of the host galaxies and the results 
of the SED fitting. 

Because the mean FWHM of the IRAC PSFs is -1.8" 
as contrasted with that of Subaru optical PSF of ~1.1", 
some objects in the IRAC images are seriously con- 
taminated by surrounding objects. For this reason, we 
checked 617 objects by eye whether they have neighbors 
in their close vicinity on the high-resolution z'-band im- 
age. Furthermore, we also examined the IRAC images 
by eye whether the neighboring objects around the LBG 
position affect the photometry to make a sample of the 
LBGs secure for the SED fitting. After these inspections, 
we selected 170 objects for subsequent analyses. In Fig- 
ure [U representative objects of the sample LBGs that 
are detected both in 3.6/im and 4.5/im bands are shown. 

The distribution of the z'-band magnitudes for the 
sample is shown in Figure [21 The z'-band magni- 
tudes for the total sample range from 26.5 to 23.5 mag 
(M X =—19.8 to —22.8 mag) with a median value of 
25.8 mag (M^^^^^=— 20.5 mag), where the absolute mag- 
nitudes are calculated with z = 4.8 and the spectral in- 
dices obtained from Iq — z' colors. In Figure [5^, the 
magnitude distribution of the sample we use in this work 
(170 objects) with that of the entire LBG sample (617 ob- 
jects) is plotted. There seems not to be much difference 
between them. We made the Kolmogorov-Smirnov test 
for the distribution of z'-band magnitude; the hypothesis 
that the IRAC LBG sample is chosen from the the same 
sample as the original LBG sample cannot be rejected at 
the 5% confidence level. We also made the test for Ic — z' 
color distribution. Again we cannot reject the hypothe- 
sis. A fraction of the uncontaminated objects in IRAC 
images in each z'-band is almost constant within la error 
and is independent of the z'-band magnitude. These tests 
support the idea that the sub-sample is not biased by our 
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Fig. 2. — (a) Distributions of the z'-band magnitudes of the sample. The di stribution of the IR AC sample used in this work is indicated 
by thick solid line. The distribution of the entire LBG sample (617 objects) bv llwata et al.l II2007I) is scaled to the total number of the IRAC 
sample (170 objects) and is indicated by dotted line. The distributions of the GOODS-N sample and the GOODS-FF are indicated by 
hatched region and cross-hatched region, respectively, (b) Distributions of z'-band magnitudes for the total IRAC sample and the samples 
categorized by the IRAC detection (from Category 1 to 4; see text for the categorization). 



selection of objects based on their lack of neighbors, i.e. 
it is randomly selected out of the total sample. Accord- 
ing to the expected fraction of the interlopers i n eac h 
z'-band magnitude bin presented bv llwata et al.l (|2007D . 
the expected number of interlopers is 25 among the IRAC 
sample of 170 objects. Among the IRAC sample, four ob- 
jects me et IRAC-selecte d extremely red object (lERO) 
criteria (|Yan et al.ll2004f ). However, even if these objects 
are truly z < 3 lEROs, they do not affect main results 
in this paper. 

The cross-matching between sources in the Subaru im- 
ages and IRAC images was made with a 1.2" search ra- 
dius. Most of the IRAC counterparts are found within 
0.2" — 0.4", which is almost comparable to the accuracy 
of our astrometry. With this search radius, 105 objects 
are detected in IRAC 3.6/im and/or 4.5/im. Among these 
objects 64 are detected both in 3.6/im and 4.5//m bands, 
and 29 and 12 are detected only in 3.6/im and 4.5/im 
band, respectively. The other 65 objects are detected 
neither in 3.6/im nor 4.5/im bands. These objects are 
grouped into four categories: objects detected both in 
3.6/im and 4.5/im (Category 1), objects detected only 
in 3.6/im (Category 2), objects detected only in 4.5/im 
(Category 3), and objects detected neither in 3.6/im nor 
4.5/im (Category 4). Note that the Category 2-4 ob- 
jects are not undetected due to the blending with the 
neighboring sources but are intrinsically faint in IRAC 
bands because we selected objects that are isolated and 
free from contaminations from nearby sources in IRAC 
bands. FigureOa shows the distributions of z'-band mag- 
nitudes of the samples categorized from 1 to 4. It is no- 
table that objects in the Category 1 sample are generally 
brighter in z'-band than those of the other categories. 
The magnitudes of Category 1 — 3 objects are presented 
in Table [Tl which will used for SED fittings. 

The distributions of the magnitudes in 3.6/im and 
4.5/im band are shown in Figure [3] and [3) In 3.6/im band, 
the magnitudes of our sample range from 7713 g^m=21.5 
to 26.9 mag, which corresponds to a range of Mr'=—2A.8 
to —19.4 mag at z = 4.8. In the 4.5/im band, the magni- 



tudes of our sample range from 7714. 5^^=21. 3 to 26.3 mag, 
which corresponds to a range of Mi/=— 25.0 to —20.0 mag 
at z = 4.8. Figure [3] shows the absence of the faint ob- 
jects in the GOODS-FF due to the shallower limiting 
magnitudes both in 3.6/im and 4.5/im bands. The figure 
also shows that the absence of the bright objects in the 
GOODS-N compared to the GOODS-FF. This may be 
due to a fact that the LBGs brighter in z'-band tend to 
reside in the GOODS-FF rather than in the GOODS-N, 
probably by chance. Figure |4] shows that the Category 2 
and 3 objects are generally fainter than the Category 1 
objects. 

3. POPULATION SYNTHESIS MODELS AND SED 
FITTING 

We then perform SED fitting by using the SEDfit soft- 
ware package^ (Sawicki, in prep.), whi ch employs es- 
sentially the same algorithm as that by ISawicki &: Yed 
(1998). 

We generated a set of model SEDs with a pop- 
ulation synthesis model as follows: We used the 
iBruzual fc CharlotI (|2003f) (hereafter BC03) population 
synthesis code wit h Padova 1994 evolutionary tracks. 
The Salpeter IMF (|Salpeteilll955h with the mass range 
of O.IMq — IOOA/q was used. Although this combination 
of the modeling might not be modern now, we intend to 
compare our results with previous studies of LBGs at 
the lower redshifts to see the evolution. Metal abun- 
dance was adopted to be 0.2^©, by considering that 
metallicity of galaxies at z = 2 — 3 are still sub-solar 
(iPettini et al. 2001; Erb et al.ll20"0 6: Halliday et al. 200i) 
and the metallicity of z ~ 5 LBGs is suggested to be at 
least Z - 0.1Z© ' ((Ando et al.|[2007[) . We adopted the 
constant star formation history (CSF). In Appendix |A| 

^ This code uses a standar d minimiza t ion al gorithm. The 
main differenc es compared to ISawicki fc Yed (119981) are that this 
code uses the ICalzetti et al.l 1)20001 ) law inst ead of the ICalzettil 
III997I) and the population synthesis model by 'Bruzual & Charlo^ 
(2003) instead of that by Bruzual & Chariot (1993). This package 
will soon be made publicly available and in the meantime can be 
obtained by contacting M. Sawicki at sawicki@ap.smu.ca. 
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Fig. 3. — Distributions of magnitudes in 3.6/im (left panel) and A.5iim (right panel) band of the sample, and those for the GOODS-N 
and the GOODS-FF region. 




3.6pm Magnitude (AB) 4.5pm Magnitude (AB) 

Fig. 4. — Distributions of magnitudes in 3.6/tm (left panel) and 4.5/im (right panel) band of the sample, and those for the Category 1, 
2, and 3. 

we investigate effects on the results by adopting other lation and is put into the model spectrum. IStark et al.l 

metallicities (I.OZq and 0.02^©) and star formation his- (|2007| ) estimated that the contribution of the Ha line to 

tories (instantaneous burst, exponentially declining mod- the flux density in the 3.6/im band is 10 — 20% for z ~ 5 

els with T—1 Myr, 10 Myr, 100 Myr, 1 Gyr, and 10 Gyr, LBGs and d oes n ot affect the stellar mass significantly, 

and two-component star formation history models). The I E vies et al.l (|2007f ) also reported the Ha contribution for 

universe at z ~ 5 is ~1.2 Gyr old and no object can z ~ 6 LBGs is < 10% and the results in their paper 

be older than that. However, we allowed for model ages do not change within errors. We examine this effect for 

of up to 20 Gyr that is the oldest age BC03 provides our sample LBGs at z ~ 5. For the z = 4.8 galaxies, 

as a check on the fits. The BC03 uses 221 age steps the effects of Ha inclusion on the magnitudes in 3.6/im- 

from 0.1 Myr to 20 Gyr, which are not equally-spaced band range from '^0.1 (for model age = 1 Gyr) to ^ 0.7 

in logarithmic scale. The SEDfit resamples this 221 to mag (for model age=10 Myr), depending on the age of 

51 equally-spaced logarithmic age steps both to speed the model spectrum. The contribution of the Ha for the 

up the calculations and to avoid having to deal with the model age of 25 Myr, which is the median value of the 

unequally-spaced scale of the original 221 models. best-fitted age as we mention in tj4.21 is ~ 0.5 mag. Since 

We took into account Ha emission line in the model this difference is larger than the typical errors in 3.6/im 

spectrum. The Ha emission line comes into the 3.6/im band and is not negligible, we take into account the Ha 

band if the redshifts of the LBGs are 4.0 < z < 5.0. The line and run the SED fitting. It should be worth empha- 

procedure of putting the Ha line in the model spectrum sizing that the inclusion of Ha emission line in the model 

is not equipped in the SEDfit. The Ha luminosity is cal- spectrum improves the fit very much without increasing 

culated from the model star forma tion rate (i.e . , intr insic the number of free parameters as shown below. Details 

star formation rate) by using the iKennicuttI (| 19981) re- for the inclusion of Ha line and effects on the results are 
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Fig. 5. — Best-fit model SEDs {open triangles) and observed SEDs {filled circles). Best-fit model spectrum is also shown with solid line. 
Although y-band magnitudes are plotted, the data were not used in the fitting procedure. In each panel, object ID and its best-fitted 
parameters are shown. 
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discussed in Appendix [Bl 

Resulting model spectrum was attenuated by internal 
dust with extinction values ranging from E{B — V) = 0.0 
to 0.8 mag in a st e p of .01 mag using the extinction law 
bv ICalzetti et al.l (]2000f). We also tested the effects of 
using the Small Magellanic Cloud (SMC), Large Mag- 
ehanic Cloud (LMC), and Milky Way (MW) extinction 
laws (see details in Appendix IXj) . Finally, the spectrum 
was attenuated by intergalact ic medium (IGM) using the 
prescription by Madau l)1995f ). 

We fixed the redshift to 2; = 4.8 in order to reduce 
the number of fr ee parameters. Acc ording to the selec- 
tion function bv llwata et al.l ()2007t ). the expected red- 
shifts of the LBGs range from z = 4.3 to z = 5.3 and 
the aver age redshift is z = 4.8. Also, the spectroscopic 
study bv lAndo et al.l ()2007[ ) shows that the distribution 
of the identified redshifts is broadly consistent with the 
expected distribution with the mean redshift of z = 4.8. 
We thus adopted the redshift of the objects in our sam- 
ple as 4.8 in the SED fitting. Note that two objects in 
the sample are spectroscopically confirmed {z = 4.70 and 
4.62), but we take the redshifts of these objects to be 4.8. 
We examine the effects of this assumption on the stellar 
mass in Appendix [Cl and find that fixing the redshift 
induces a systematic offset of only 0.06 dex and the un- 
certainty is 0.22 dex compared to making the redshift 
free. We additionally assess the uncertainty by assigning 
the redshift r andom ly along the expected distribution by 
llwata et al.l ()2007l ). The uncertainties are 3 — 5 times 
smaller than the errors in the SED fitting (see Appendix 
p. 

In this SED fitting analysis, except for a testing case 
of redshift free fitting, we used the data in Ic, z' , IRAC 
3.6/im, and 4.5/j,m bands. We did not use F-band mag- 
nitude to avoid the uncertainty due to the fluctuation 
of IGM absorption. For Category 2 and 3 objects, the 
upper limit magnitudes were not used in the SED fit- 
ting. The free parameters in the SED fitting are age, 
color excess, and scaling normalization (stellar mass and 
SFR). 

4. RESULTS 

4.1. Stellar Mass 

We select via eye inspection 170 objects that are iso- 
lated and are not contaminated by neighboring objects 
in 2:'-band and IRAC images out of the entire LEG sam- 
ple consisting of 617 objects. We made the SED fitting 
for 105 objects that are detected in IRAC 3.6/Ltm and/or 
4.5/im among the 170 objects. We did not make the 
SED fitting for the other 65 objects detected neither in 
3.6/im nor 4.5/im. Figure [5] shows representative exam- 
ples of the observed SEDs and the best-fitted models 
for 12 objects (6 objects from the Category 1, 3 objects 
from the Category 2, and 3 objects from the Category 
3). The best- fitted parameters (stellar mass, stellar age, 
color excess, and star formation rate) of each object are 
summarized in Table [21 Almost all of the SEDs are well 
reproduced. The effect of including Ha emission can be 
seen particularly in the second row of Figure [5l the ex- 
cess due to the Ha emission is significant in the IRAC 
3.6/im band. However as seen in the bottom row of Fig- 
ure [U for the Category 3 the model flux density is larger 
than the upper limit in the 3.6/im band. This may be 



caused by the assumption of z = 4.8. Category 3 objects 
may show an excess in IRAC 4.5/im band rather than 
in 3.6/im. Since the Ha emission comes into the 4.5/im 
band if z > 5.1, the redshifts of these LBGs may be larger 
than 5.1. In fact, if we take the redshift as a free param- 
eter in the range from z = 3.8 to 5.5 with Az = 0.1, 
the value of the chi-square is minimum a,t z ^ 5.2 or 
is almost comparable to those at z < 5.0, suggesting the 
presence of Ha emission in the 4.5/i,m band. We proceed, 
however, assuming that all of our sample objects are at 
z = 4.8. 

The distribution of the best-fitted stellar masses is 
shown in Figure [6^. The derived stellar masses of the 
whole sample galaxies range from 10^ to IO^^Mq with a 
median value of 4.1 x 1O^M0. The typical error at the 
90% confidence level in each SED fitting is ~0.4 dex. As 
we discuss in Appendix \^ the differing star formation 
histories, metallicities, and dust extinction laws affect 
the output parameters of the fitting to some degree. Age 
and color excess are most affected and stellar mass is 
least affected. The uncertainties of these effects on the 
stellar mass are ~0.6 dex at most. Figure [6^ indicates 
that some massive galaxies (> IO^^Mq) have already 
been assembled at z ^ 5 when the universe was only 
1.2 billion years old. The best-fitted stellar ages of 20 
objects are, however, older than the cosmic age (~1.2 
Gyr) at z '--^ 5. In the distributions of Figure [6l these 
objects are shown as shaded regions. It is noteworthy 
that the stellar masses of these overage objects are typi- 
cally large. A cause for the large ages is considered to be 
due to the assumption of constant sta.r form ation history 
(jSawicki fc Yeelll998t iPapovich et al.ll200lD . In general, 
the derived ages are older than those by assuming other 
star formation histories such as an instantaneous burst 
or an exponentially declining star-formation history. Al- 
though we focus on results for objects whose best-fitted 
ages are younger than 1.2 Gyr, the presence of this pop- 
ulation should be kept in mind. We will come back to 
this problem when we derive the stellar mass function 
and density. 

A relationship between magnitudes in the 4.5/im band 
(which corresponds to approximately the rest-frame i' 
band) and stellar masses is shown in Figure [71 A clear 
correlation between the 4.5/im magnitudes and the stel- 
lar masses (log(M*/M0) = -0.56 x m4.5^,„ + 23.27) 
is seen (especially for Category 1), indicating that the 
rest-frame optical flux is a good indicator of the stellar 
masses. The correlation is not exactly linear; the mass- 
to-light ratio is larger in the brighter objects. 

The FigurelHi shows that the objects in the GOODS-N 
are relatively less massive than those in the GOODS-FF. 
The median stellar masses of the objects in the GOODS- 
N and the GOODS-FF are 2Ax10^Mq and l.OxlOi°M0, 
respectively. It is reasonable because in the GOODS- 
FF the LBGs faint in IRAC bands are absent, while the 
bright ones tend to reside in. The number of massive (> 
10^°Mq) objects in the GOODS-N is smaller than that 
in the GOODS-FF. This deflcit of the massive objects in 
the GOODS-N may be due to cosmic variance. Figure 
[8] shows distributions of stellar masses for the Category 
1, 2, and 3 objects. The median stellar masses of the 
Category 1, 2, and 3 objects are 6.9 x IO^Mq, 2.9 x 
IO^Mq, and 3.6 x IO^Mq, respectively. As a whole, the 
Category 2 or 3 objects are relatively less massive than 
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the Category 1 objects. This is also reasonable if we 
recall the faintness of the Category 2 and 3 objects in 
the rest-frame optical wavelength (Figure . 

As shown in Figure [21 there seems to be no clear cor- 
relation between the rest-frame UV absolute magnitudes 
(uncorrected for dust extinction) and the stellar masses 
in LBGs at z ~ 5. However, we compute a correlation 
coefficient for the relation, and obtain r = —0.3. We 
can reject the null hypothesis that there is no correla- 
tion between the UV absolute magnitude and the stel- 
lar mass at the 5% confidence level. The median stel- 
lar masses in 0.5 mag bins show a loose correlation be- 
tween the rest-frame UV absolute magnitudes and the 
stefiar masses: log(M*/MQ) = -0.38 x M^^^^^ + 1.64. 

IShaplev et al.l (|200lL l2005( l found no correlation between 
the UV absolute magnitude and the stel l ar m ass for 
LBGs at z = 2 - 3. However. fSawicki et "all (|2007t ) found 
that sub-L* LBG s at z ~ 2 tha t are m uch fainter than 
those st udied bv IShaplev et al.l ()2005f ) show the corre- 
lation. iPapovich et al.l i 2001h also found the correla- 
tion between the UV absolute magnitude and the stellar 
mass for LBGs at z ~ 3 by using a deeper sample than 



that by IShaplev et al.l (|2001h . The correlation between 
the rest-frame UV magnitude and the stellar mass for 
z ^ 5 sample is mu ch weaker than that found at z ~ 2 
(jSawicki et al.ll2007t Sawicki, in prep.). 

4.2. Ages, Dust, and Star Formation Rates 

The distributions of the best-fitted ages and color ex- 
cesses are also presented in Figure and [nt, respec- 
tively. The typical error of the age for each object is 
~1.0 dex. The median value of the ages estimated for 
our sample is 25 Myr. The median values of the ages es- 
timated for the GOODS-N and the GOODS-FF sample 
are 19 Myr and 31 Myr, respectively. The median values 
of the ages of the Category 1, 2, and 3 objects are 35 
Myr, 25 Myr, and 5 Myr, respectively. The typical error 
of the color excess for each object is '^O.l mag. The me- 
dian value of the color excesses estimated for the total 
LBG sample is 0.22 mag. The median values of the color 
excesses estimated for the GOODS-N and the GOODS- 
FF sample are 0.18 mag and 0.25 mag, respectively. The 
median values of the color excesses of the Category 1, 
2, and 3 objects are 0.20 mag, 0.20 mag, and 0.33 mag, 
respectively. 
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4.5Mm Magnitude (AB) 

Fig. 7. — 4.5/jm magnitudes vs. the derived stellar masses. Cat- 
egory 1 objects and Category 3 objects are indicated by circles and 
squares, respectively. The vertical error bars are 90% confidence 
level and the horizontal bars show la errors in magnitudes. 



It is known, however, that rest-frame UV-optical colors 
are degenerate with respect to age and dust extinction. 
Whether the colors are explained by extinction or age 
is hard to be specified uniquely. Generally speaking, for 
LBGs at z ~ 5, this age-dust degeneracy may be broken 
by adding NIR data to the SEDs, and we may be able to 
improve our estimation of a ges and dust conten t. To test 
this, by using the sample bv 'Stark et al.l (|2007f ) including 
objects with spec-z and phot-z, for which both J and Ks 
data are available, we compare results derived by using 
the J and Ks data and those without using the data. 
The details are presented in Appendix [D] The results 
show that differences with and without NIR data are not 
large, although error bars are large. The stellar masses 
with and without the J and Kg data agree with each 
other within a factor of ~ 3. However, the errors in the 
J and Ks data used in the test are generally larger than 
those in other bands and the weights of the NIR data to 
the SED fitting are relatively small. Thus this may cause 
the small differences in the test. Therefore, sufficiently 
deep NIR data are desirable to better constrain age and 
color excess. 

The distribution of the best-fitted SFRs is presented in 
Figure [6]i. The typical error of the SFR for each object 
is ~0.5 dex. The median value of the SFRs estimated 
for the total LBG sample is 141M0yr~^. The median 
values of the SFRs estimated for the GOODS-N and the 
GOODS-FF sample are lO4M0yr-i and 191Moyr-\ re- 
spectively. The median SFRs of the Category 1,2, and 3 
objects are UOMqyi-^, lllMoyr^i, and 1023Moyr~i, 
respectively. Figure [6]i shows the existence of galaxies 
which show high SFRs. This is because of the existence of 
a large amount of dust as presented above. T he apparent 



SFR derived from L* at z ~ 5 derived bv llwata et al.l 

(12001 using lMadaul ^ (|1996[ ) relat ion is ^ 2OM0yr ^. If 
we use the median value of color excess of the galaxies of 
0.22 mag, the extinc t ion co rrected SFR is ~ 16OM0yr~^ 
using [Calzetti et"aLl (|2000f ) extinction law, which is con- 



sistent with the value derived from SED fitting. 

5. DISCUSSION 

In §31 we showed the stellar properties of LBGs at 
z ~ 5 derived from SED fitting. In this section, we com- 
pare the results in this work with previous studies for 
galaxies at other redshifts. We also construct a stellar 
mass function of LBGs at z ^ 5 and derive the stellar 
mass density at z ~ 5. 

5.1. Comparisons with LBGs a< z = 2 — 3 

First we compare the properties of LBGs at z ~ 5 
with those at z = 2 — 3 (cosmic age of 3.2 Gyr— 2.1 
Gyr). Her e we use the term inology of LBGs including 
BM/BX ( Stei del et aI]|2004D . The distributions of the 
output parameters for our sample are compared with 
those of z = 2 — 3 samples in Figure [TUl where the 
histograms are normalized so that its peak value equals 
unity f or comparison. For sa mple LBGs at z = 2 and 3 , 
we use IShaplev et ahl ([2005') and IShapleveTall ((200lfl . 
respec tively. All three sample s are fitted using mod- 
els bv lBruzual fc Chariot! (|2003f ) with the Salpeter IMF 
consta nt star formation history, and the iCalzetti et al.l 
()2000( ) extinction law. 

In order to compare the samples fairly, their faintest 
UV luminosities at 1500A are on an equal footing with 
those of ours. The distributions of the rest-frame UV 
absolu te ma gnitudes o f our z ~ 5 sample, the z = 3 sam- 
ple bv | Shap lev et al] (|2001l ). and the z = 2 sample by 
IShaplev et al.l ^ (^2005) are presented in Figure [TT] While 
the UV absolute magnitudes of the z = 3 and the z = 2 
samples lie on the range of —19.3 mag to —21.7 mag and 

— 19.6 mag to —22.5 mag, respectively, that of our sample 
ranges from —19.7 mag to —23.0 mag. We use the sam- 
ple galaxies whose UV absolute magnitudes are brighter 
than —19.7 mag. Figure [121 shows the distribution of the 
rest-frame optical (5500A) a bsolute magnitudes o f the 
z ^ 5 sample, z = 3 sample b ylShapl ev et all (|200lD . and 
z = 2 sample bv lShaplev et a l. (2005). While the faintest 
magnitudes of these samples are almost the same (from 

— 19 to —20 mag), the brightest magnitudes are some- 
what different. The optical absolute magnitudes of the 
z = 5 sample range from —19.1 mag to —24.9 mag, while 
those of the z — 3 and z = 2 samples range from —19.7 
mag to —22.8 mag and from —20.3 mag to —23.8 mag, 
respectively. There are clear deficits in the bright parts 



of both Muv and M, 



^optical 



distributions of the z = 3 sam- 



ple compared to the z ~ 5 sample. This is probably due 
to a smaller survey volume for the z = 3 sample than 
that for the z ^ 5 sample. The differences may also be 
attributed to cosmic variance. 

In the upper right panel of Figure [TTl the stellar masses 
against the rest-frame UV absolute magnitudes of the 
z = 2, z = 3 samples and our z 5 sample are plotted. 
The relation between the stellar mass and the UV abso- 
lute magnitude varies from Z'^5toz = 2 — 3 toward 
large masses at a fixed UV absolute magnitude, although 
the correlation is not strong (see details in M.ip . In the 
upper right panel of Figure [121 the stellar masses against 
the rest-frame optical absolute magnitudes of the z = 2, 
z = 3 samples and our z ^ 5 sample are plotted. The 
correlation between the stellar mass and the optical ab- 
solute magnitude also varies toward large masses at a 
fixed UV absolute magnitude from z~5toz = 2 — 3. 
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Fig. 8. — Distributions of the best-fitted parameters. The distributions of the best-fitted stellar masses, ages, color excesses, SFRs are 
plotted in panel (a), (b), (c), and (d), respectively. The LBG sample is divided into four categories (see text). Shaded regions indicate the 
overage LBGs; these overage histograms are plotted cumulatively on top of the < 1.2 Gyr histograms. See text for details. 



In Figure [TUk . the distribution of stellar masses for 
LBGs at z '--^ 5 superposed on those at z = 2 and 3 
is presented (The comparison of the distribution of stel- 
lar masses is also presented in the upper left panels of 
Figure [TT] and Figure fT2|) . Figure [TOk shows that the 
stellar masses of LBGs at z ~ 5 are smaller than those 
at z = 2 — 3 on average. While the median of the stellar 
masses of the z ~ 5 LBGs is 4.1 x lO^M©, the medians 
of the masses of the z = 2 and 3 LBGs are 1.7 x 10^°Mq 
and 1.3 x lO^^Af©, respectively. Therefore, the median 
stellar mass of z 5 LBGs are smaller by a factor of 
3-4 than that of z = 2-3 LBGs. Note t hat while 
iShaplev et all (|200lD and iShaplev et all (|2005l ) used so- 
lar metallicity models {I.QZq), we use sub-solar metal- 
licity models (O.2Z0). If we assume the metallicity of 
1.0.^0, the stellar mass decreases by a factor of 1.2, and 
the diffe rence between the distributions is even more sig- 
nificant. IVerma et al.l ([2007) also found that the typical 
stellar mass of z ^ 5 LBGs is 5 — 10 times lower than the 
z = 3 LBGs. The stellar masses of these z ~ 5 LBGs are 
almost compar able to those of muc h fainter (suh-L^y) 
LBGs at z 2 (jSawicki et al.ll2007t Sawicki, in prep.). 



As a whole, star formation ages of z ^ 5 LBGs are 
younger than those of z = 2 - 3 LBGs (Fig. ]Wjp). While 
the median age of our sample is 25 Myr, the median ages 
of z = 2 and z = 3 LBGs are - 600 Myr and - 300 Myr, 
respectively. The bimodal distributi on in z = 2 — 3 sam - 
ples is not seen in the z 5 sample. IVerma et al.l (|2007l l 
also found that z ~ 5 LBGs are significantly younger 
than the z — 3 LBGs. However, typical ages are < 10 
Myr, which are about one order of magnitude smaller 
than our results. The cause of this difference is not clear. 
The distribution of color excess of the z ~ 5 LBGs seems 
to suggest that the amount of dust extinction in z 5 
LBGs may be slightly larger than that of the z '--^ 3 
LBGs and similar to that of the z 2 LBGs (Figure 
\Wb). The median color excess of our sample is 0.22 mag, 
and the median color excesses for z = 2 and 3 LBGs 
are 0. 20 mag and 0.16 mag, respectively. IVerma et al.l 
(|2007f ) also found that z ~5 LBGs have a typical color 
excess of ^ 0.2 mag. The SFR is higher than those of 
z = 2 - 3 LBGs (Fig. [TOH). While the median SFR 
of z 5 LBGs is 141MQyr^^, the medians of SFRs are 
52M0yr^^ and 43MQyr^^ for z = 2 and z — i samples, 
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Fig. 9. — Absolute UV magnitudes (uncorrected for dust extinc- 
tion) vs. stellar masses for Category 1 objects (circles), Category 2 
objects (triangles), and Category 3 objects (squares). Large open 
circles show median stellar masses in 0.5 mag bin. Solid line is the 
regression line for these points. 



respectively. Specific SFR is al so larger in z ~5 LBGs 
than that in z = 2 — 3 LBGs. iVerma et all (|2Q07l ) also 
found that the typical SFR (~ 500 Moyr"!) of z - 5 
LBGs is ^ 10 times higher than the z = 3 sample. From 
these comparisons, we suggest that galaxies at z 5 are 
forming stars very actively, and in consequence, they are 
dusty and we may see the early phase of these activities; 
we may witness the evolution of stellar populations of 
galaxies from z = 5 to z = 2. 

In the bottom sub-panel of Figure [TUk . we show dis- 
tribution of stellar masses of z = 3 LBGs assuming that 
each galaxy of our sample keeps the SFR derived from 
the SED fitting until z = 3. The distribution shifts to- 
ward larger mass and the median value of the distribu- 
tion is 1.3 X 10^"'^ M0. Likewise, we plot the expected 
distribution of stellar masses of z = 2 LBGs in the top 
sub-panel of FigurefTUk. Again, the distribution shifts to- 
ward larger mass than observed at z = 2, and the median 
value of the distribution is 2.9 x 1O^^M0. Thus, these 
suggest that the SFR decreases from z '--^ 5 to 3 — 2, pro- 
vided that z = 2 — 3 LBGs are direct descendant of the 
z ~ 5 LBGs. Alternatively, the descendant of the z ~ 5 



LBGs may be massive objects at z 
or sBzKs. 



2 — 3 such as DRGs 



5.2. Comparisons with LBGs at z = 5 — 6 

We compare the properties of our sample galaxies with 
those at z = 5 — 6. For th e z = 5 samples, we us e 
samples by IStark et all (I2007D an d IVerma efall (|2007D . 
The sample bv IStark et al. (|2007t ) consists of 14 spectro- 
scopically confirme d objects at z 5 and the sample by 
IVerma et al.l (|2007[ ) consists of 21 ^-dropouts, 6 of which 
are confirmed to be at z 5 by spectroscopy. For the 
z = 6 samp l es, we use samples by ma,n et al. (200^ and 
lEvles et all (I2007D. The s ample bv lYan et all (|2006f) and 
that bv lEvles et al] (|2007f ) consist of 53 and 9 i'-dropouts, 
respectively. 



For fair comparison, we checked the limiting magni- 
tudes for these samples both in the rest-frame UV and 
in the rest-frame optical band. The rest-frame UV and 
optical magnitude ranges of the other studies are almost 
the same as ours. The faintest limits of the absolute 
magnitude in the rest-frame UV and the rest-frame op- 
tical are M s ~ —20 mag and M s ~ —20 mag, 

ISOOA ° 5500A °' 

respectively. 

In Figure [T31 we plot the distributions of derived pa- 
rameters for the various samples, where the peaks of the 
histograms are normalized to be unity for comparison. 
Note that because the age and the color excess are not 
shown explicitly in the paper bv lYan et al.l (|2006[ ). we do 
not plot them in the figure. 

As illustrated in Figure [T3k . the range of the stellar 
masses for our sample agrees with other studies. The 
stellar masses are widely distributed from 10® M© to 
10^^ Mq. While the median stellar m ass of our sample is 
4.1 X lO^M© those o f the sample s of iSj arkjF al. (200?;) 
IVerma et all (|2007[ ). lEvles et al.l (|2007D . and lYanet alj 
(1200^ are 7.9 x lO^M©, 2.0 x IO^Mq, 1.6 x 10^°Mq, 
and 9.6 x lO^M©, respectively. It seems to be somewhat 
strange that the representative stellar mass for the z = 5 
LBGs is less massive than that for the z = 6 LBGs. It is 
also noteworthy that the models used in the SED fitting 
for the z = 6 sample are slightly different from those we 
used; we will discuss the effects of differing models below. 

Figure [T5b shows the comparisons of the distributions 
of the star formation ages for our sample and other stud- 
ies. The resulting age distribution of the other studies 
tend to be younger or older than our sample: While the 
median value of the age of our sample is 2 5 Myr, and 
the median ages for the z = 5 sa mple by IStark et al.l 
(|2007t l and for the z = 6 sample by lEvles et all (|2007l ) 
are 255 Myr and 400 Myr, respectively. On the other 
hand, the ages for the z = 5 sample by IVerma et al.l 
()2007t l are typically younger than 10 Myr. In Figure IT3h. 
we plot the distribution of color excesses for our sample 
and t he other studies. The derived color excesses for 
z = 5 (jStark et al.l l2007^ and z = 6 samples (jEvles et al.l 
|2007| ). most of which are close to zero (the median val- 
ues are 0.01 mag a nd 0.00 mag for lStark et al] ()2007D and 
lEvles et al] ()2007( ). respectively), are much smaller than 
our result of the median value of 0.22 mag. In contrast, 
the resulting color excesses for our sampl e are almost 
compa rable or slightly smaller than those bv lVerma et al.l 
(|2007h . 

In Figure [T3H. we plot the distribution of SFRs for our 
sample and the other studies. The SFRs for our sample 
are much higher t han those for z = 6 LBGs ()Yan et al.l 
120061 lEvles et al.l |2007.) but slightly lower than those 
of the z = 5 sample by IVerma et al.l ()2007l ). The me- 
dian values of SFR s for our sample and the sample by 
IVerma et all (|2007l ) are 141M0yr-i and 4OOM0yr-\ re- 
spectively. On the other hand, the median SFR of the 
z = 6 samples is < lOMQyr^^. This difference between 
z = 5 and z = 6 samples may be due to the difference 
of color excesses. Since extinction in z = 6 sample is 
negligible, the extinction corrected SFR is low. Mean- 
while, for the z = 5 samples, the moderate amount of 
dust extinction makes the intrinsic SFR higher. 

It is worth noting that model ingredients in the SED 
fitting of these studies for deriving the parameters are 
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Fig. 10. — Distri butions of the best-fitted parameters of our 2 ~ 5 sample ( thick line) and those of z = 2 sample (hatched) from 
IShaplev et al.l I I2005I ) and z = 3 sample (cross-hatched) from lShaplev et al.l I I2001I) . The distributions of the best-fitted stellar mass, age, 
color excess, and star formation rate are plotted in panel (a), (b), (c), and (d), respectively. In panel (a), dotted lines of top and bottom 
sub-panels indicate distributions of expected stellar masses at z = 2 and 2 = 3, respectively, if our sample galaxies at z ~ 5 keep the SFRs 
derived from the SED fitting. For comparison, the peaks of all distributions are normalized to unity. 



different. Thus, it is necessary to take into account these 
effe cts when we comp are o ur results with ot her stud- 
ies. IStark et all (|2007D and lEvles et al.l (|2007D used the 
BC03 with th e Salp eter IMF and the extinction law by 
iCalzetti et al.l pOOu) , but they assumed solar metallicity 
(I.OZq) and the various star formation histo ries. We re- 
fit the observed SEDs of the z = 5 objects bv lStark et al.l 
()2007() . including 14 objects with sp ec-z and 59 o bjects 
with phot-z, and z = 6 objects by lEvles et all ()2007D 
with the same models as we used for our sample: BC03, 
Salpeter IMF, 0.2 Z(t,, constant s t ar for mation, and the 
extinction law by ICalzetti et al.l ()2000[ ) . The resulting 
distributions are presented in Figur e The refitted 
parameters for the z = 5 sample by IStark et al.l ()2007f ) 
are somewhat different from the original re sults, while 
the re fitted results for the z = 6 sample by lEvles et al.l 
(l2007l) are almost the same as the original results. The 
median stellar masses of the the z = 5 sample and the 
z = 6 sample are 2.3 x IO^Mq and 1.8 x IO^^Mq, re- 
spectively. They decrease by a factor of 3.4 and increase 



by a factor of 1.1 from original results, respectively. The 
star formation ages also vary from original results. The 
refitted ages for the z = 5 sample decreases from the 
original result by a factor of ~10 (the median value of 
25 Myr) and is comparable to our result. On the other 
hand, the age for the z — 6 sample increases by a factor 
of 1.7 from the original result (the median value of 700 
Myr). Note that when the constant star formation his- 
tory is assumed, ages of some objects in the z = 6 sample 
exceed the cosmic age at z = 6 (~ 0.9 Gyr). While the 
color excesses of both the z = 5 and z = 6 samples are 
~ mag, the refitted color excess of the z = 5 sample is 
0.17 mag, which is comparable to our results, and that 
of the z = 6 sample is 0.0 mag. The median values of 
the refitted SFR are SGMgyr"! and 2OM0yr-i for the 
z = 5 and z = 6 sample, respectively, and are lower than 
our result. 

Even though the model ingredients in the SED fitting 
are the same, the resulting stellar mass, age, color excess, 
and star formation rate of the z = 6 samples are different 
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Fig. 11. — Upper right: Rest-frame UV (1500A) absolute magnitu des (u ncorrected for dust extinction) vs. stellar^ masses for 2 ~ 5 
sample (filled circles), 2 = 3 sample (open triangles) by|Shaplcy et al] 11200 11 1, and 2 = 2 sample (open squares) by Shaplov et al. (2005). 
Lower right: Distributions of the rest-frame UV absolute magnitudes of the 2 ~ 5 sample {thick line), z = 3 sample (cross-hatched) , and 
2 = 2 sample (hatched). Upper left: Distributions of the stellar masses of the 2 ~ 5 sample (thick line), 2 = 3 sample (cross-hatched), and 
2 = 2 sample (hatched). 



from those of our z — 5 sample. This may imply that 
there is a significant evolution of steUar population in 
galaxies from z — 6 to z — 5. Although the time interval 
between z = 6 and z — 4.8 is just ~ 0.3 Gyr, the galaxies 
may evol ve drastically aft e r the end of reioniz ation epoch 
at ^ - 6 (|Fan et all 120061 : 1 Tot ani et al.ll2006f l. However, 
the small sample size of the previous studies prevents 
us from concluding that the differences of the resulting 
parameters in the SED fitting between z = 5 and z ~ 6 
samples are significant. 



5.3. The Stellar Mass Function of LBGs at z ^ b 

The large sample of LBGs whose stellar masses are es- 
timated robustly allows us to derive the stellar mass func- 
tion of LBGs at z ~ 5. Our sample is originally selected 
in the optical band, and is selected as uncontaminated 
objects in IRAC images by eye inspection. The sample is 
also affected by incompleteness in IRAC images. There- 
fore, we estimate the number density per log(A/*/MQ) 



as follows: 
0(/o5(Af,/Af0)) 



t,3 



eifdetv-ffAlog{MjMQ) ' 



(2) 

where i and j are bin numbers for z'-band and 4.5/im- 
band magnitudes, respectively, and Nij is the number 
of objects entering each z'-band bin and 4.5/im-band bin 
in the log{M^,/MQ) bin. /™* refers to a fraction of in- 
terlop ers estimated in i-th z'-band bin by Iwat a et al.l 
()2007f ). Z'"^' is a fraction of uncontaminated objects in 
IRAC images. As we discussed in H2.31 the percentage 
is independent from the z'-band magnitude, fp*' is a 
detection rate in the IRAC 4.5^m-band. For the Cate- 
gory 3 objects, which are undetected in 4.5/im-band, we 
use a detection rate in 3.6/xm-band instead of that in 
4.5/im-band. V^^^ is an effecti ve volume as a func tion of 
z'-band magnitude taken from llwata et al.l (|2007t l. 

Wc found that 20 out of 105 objects are best-fitted 
with models which have larger age than the cosmic age 
at z ~ 5. It is important to see the contribution from 
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Fig. 12. — Upper right: Rest-frame optical (5500A) absolute magnitudes (uncorrected for dust extinction) vs. stellar masses for z ~ 5 
sample (filled circles), z = 3 sample (open triangles) by Shaplcv et al. (200l|), and z = 2 sample (open squares) by Shaplcy et al. (2005). 
Lower right: Distributions of the rest-frame optical absolute magnitudes of the z ~ 5 sample (thick line), 2 = 3 sample (cross-hatched), 
and 2 = 2 sample (hatched). Upper left: Distributions of the stellar masses of the 2 ~ 5 sample (thick line), 2 = 3 sample (cross-hatched), 
and 2 = 2 sample (hatched). 



these objects to the steUar mass function, because these 
objects have generally large stellar masses, and their in- 
clusion would affect the massive end of the stellar mass 
function. In Figure [151 we plot the stellar mass func- 
tion by gathering the Category 1, 2, and 3 objects. The 
stellar mass function excluding the overage objects is in- 
dicated by filled circles, and that including the objects is 
indicated by open circles. The error bars of the number 
densities are Poisson errors^. While the number densities 
including the overage objects are higher than those with- 
out using the overage objects by a factor of 3 — 4 in the 
massive part lO^M© - lO^-^M©), the effect of the 

''' We estimated the error arising from the uncertainty of the stel- 
lar mass derived in the SED fitting by the re-sampling method. The 
stellar mass function is re-calculated by using the stellar masses 
that perturbed randomly by Gaussian distribution with Ic error 
of the SED fitting. Repeating this process 10000 times, we esti- 
mate la error of the number density in each mass bin. The average 
error is 0.17 dex, which is 1.5 times larger than the average Pois- 
son error, but is generally smaller than the uncertainties from the 
choice of the star formation history. 



overage objects is small in less massive part (< 10^^ Mq). 
If we restrict the maximum age of the models to 1.2 Gyr 
in the SED fitting, the main results, especially of stel- 
lar mass, do not change so much. In this case, the stellar 
masses are affected only in the massive part. In the most 
massive bin of the stellar mass function, the number den- 
sity decrease by a factor of ^ 3. In the other mass bins, 
the changes are within a factor of 1.4. 

As we discuss in Appendix]^ the most uncertain factor 
in deriving the stellar mass is the choice of the star forma- 
tion history. We examined the instantaneous burst, ex- 
ponentially declining history, and two-component models 
as well as constant star formation history. Wc found that 
there is no systematic offset between stellar mass derived 
with different star formation histories, but there are scat- 
ters of cr = 0.6 and 0.3 dex in the case of instantaneous 
burst and exponentially declining models, respectively, 
and also we found that, in the case of two-component 
models, the stellar masses of some objects increase by 
~ 1 dex but those of the majorities are comparable to 
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Fig. 13. — Distributio ns of the best-fitted parameters of our z ~ 5 sample and those of z = 5 — 6 samples from the Uterature ("Sta rk et al.l 
II2007I ) and lVerma et al . (2007) for z ~ 5, and lYan et al.l I I2006I) and lEvles et al.l I I2007I) for 2 = 6). The distributions of the stellar mass, age, 
color excess, and star formation rate are plotted in panel (a), (b), (c), and (d), respectively. For comparison, the peaks of all distributions 
are normalized to unity. 



those in the case of constant star formation history. In 
Figure [TSl we plot the stellar mass functions by assum- 
ing the constant star formation, instantaneous burst, and 
exponentially declining histories. Although it seems that 
the choice of the star formation history does not affect 
the stellar mass function as a whole, in the most massive 
part, the number density adopting instantaneous burst 
is larger than that adopting other star formation histo- 
ries by a factor of 2 — 3. The stellar mass function in 
the least massive part (M* < IO^Mq) decreases regard- 
less of the adopting star formation history. Although we 
apply the completeness correction, this decrease in the 
less massive end is probably due to the limitation of our 
original sample of z' < 26.5 mag as described below. 

We discussed the stellar mass function for the sample 
galaxies which are detected in 3.6/im and/or 4.5/Lim. As 
we mentioned in about 40% of our sample of 170 

objects are detected neither in 3.6^m nor 4.5^m. These 
Category 4 objects are thought to be less massive than 
IRAC detected objects (Category 1, 2, and 3). However, 
since the Category 4 objects are not detected both in 



3.6/im and 4.5/im bands, we can not derive the stellar 
masses for these objects but can only constrain the upper 
limit of the stellar mass. If we assume the correlation 
between 4.5//m magnitude and stellar mass in Figure [7l 
the upper limits (3(t) on the stellar mass for objects in the 
GOODS-N and GOODS-FF are 1.6 x IO^Mq and 8.9 x 
IO^Mq, respectively. In Figure[T5l we show the expected 
number densities including the Category 4 objects by red 
horizontal bars. Here we assume that all of the Category 
4 objects have these limiting stellar masses. We should 
emphasize that these Category 4 objects only affect the 
less massive part of the stellar mass function. 

Our sample objects are selected with a criterion of 
z' < 26.5 mag, thus objects with z' > 26.5 mag are 
missed. We roughly estimate the contribution from the 
objects that are faint in z'-band. As we discussed in ij4.11 
the rest-frame UV absolute magnitudes are roughly cor- 
related with the median stellar masses. By using this 
relationship and t he UV luminosi ty function (UVLF) of 
LBGs at z '-^ 5 (Iwa ta et al.ll2007[ ). we estimate the stel- 
lar mass function. The resulting stellar mass function is 
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Fig . 14. — Distrib utions of th e best- fitted parameters of our z ~ 5 sample and those of z = 5 — 6 samples from literature fStark ct aU 
II2007I ) for 2 = 5 and lEvles et al.l (|2003) for 2 = 6). The best-fitted parameters of samples bv lStark et al.l (|2007l ) and lEvles et al.. (.2007. ) are 
derived by fitting with the same SED models as we used in this paper. The peaks of all distributions are normalized to unity. 



plotted in FigurefTSlas a solid curve. It disagrees with the 
observed stellar mass functions, especially in the massive 
part, because the scatter of the correlation is large and 
asymmetric; while there are some objects with faint UV 
magnitudes and large stellar masses, there is no object 
with the bright UV magnitude and small stellar mass as 
seen in Figure O Nevertheless, this hints that the fainter 
LBGs {z' > 26.5 mag) contribute to the less massive part 
(Af, < 10^ Mq) of the stellar mass function. 

Here we adopt the stellar mass function derived with 
constant star formation models and with the objects 
whose age is younger than 1.2 Gyr as a fiducial stel- 
lar mass function. The fiducial stellar mass function is 
forced to be fitted by Schechter (1976) function. We ex- 
clude mass bins of log(M,/MQ) = 8.0 and 8.5 from the fit 
because we probably largely underestimate the number 
densities at the mass bins as we mentioned above. The 
best fitted parameters^ are log(M;/M0) = 13.81^^;?^, 
0* = O.eOtoil X 10-'^Mpc"Vlog(M,) , and a* = 

* The uncertainties are 68% confidence level. 



— 1.83j;Q itg. The stellar mass function from our sample 
is fitted by only power-law component of the Schechter 
function in the mass range of log(M*/Mo) = 9.0 to 11.5. 

The derived stellar mass function for our sample is 
compared with other observations. In Figure 1161 the 
fiducial stellar mass function from our sa mple and the 
result s from photo-z selected samp les by iDrorv et al.l 
(|2005D are plotted. The samples by Drorv et al.l (|2005f ) 
are an /-selected sample (/ < 26.8 mag) in the FORS 
Deep Field (FDF) and a i^j-selected sample {Kg < 25.4 
mag) in the Great Observatories Origins Deep Survey- 
So urth (GOOpS-S). T he stellar masses of the samples 
bv IDrorv et all (|2005l) are derived without IRAC data. 
As illustrated in Figure [TBI the stell ar mass func t ion o f 
our sample agrees with the results bv lDrorv et al.l ()2005f ) 
in the most massive end. However, in most of the mass 
range, the number densities of our sampl e are signifi- 
cantly smaller than those from the sample bv lDrorv et al.l 
(2005), even though we take into account the contri- 
bution of the Category 4 objects. Our stellar mass 
function is also compared with a result by I Eisner et al.l 
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Fig. 15. — Stellar mass functions of LBGs at 5; ~ 5 (totaling cate- 
gories 1-3 objects). The stellar mass functions derived by assuming 
constant star formation (circles), instantaneous burst (triangles), 
and exponentially declining (squares) models are plotted. The stel- 
lar mass functions excluding objects whose ages are fitted to be 
older than the cosmic age at 2 = 5 are indicated by filled symbols 
and the functions including these overage objects are indicated by 
open symbols. For clarity, the data points (filled-squares, filled- 
triangles, open-circles, and open-squares) are shifted horizontally 
by -1-0.1, +0.05, —0.05, and —0.1 dex, respectively. The expected 
number densities of the sample including the Category 4 objects are 
indicated by horizontal bars (see text for the detailed derivation). 
The solid curve represents the stellar mass function derived by as- 
suming UV luminosity function of z ^ 5 LBGs and UV luminosity 
and stellar mass relation (see section [5. 3 II . 




Fig. 16. — Stellar mass function of our LEG sample at 2 --^ 5 
(circles) and that of a s ample by Eisner ct al. (2008) (squares), 
that of an FDF sample bv lDrorv et al. (2005) (filled triangles), and 
that of a GOODS sample by Drorv ct al. (2005) (open triangles). 
For cl arity, the data points by Drory et al. (2005) and Eisner ct al. 
1I2OO8I) are shifted horizontally by -fO.l and —0.1 dex, respectively. 
The expected stellar mass function from the UVLF is indicated by 
a solid line (see text for the detailed deri vation) . The theor etical 
predictions of semi-analytic al models by IBower et al.l p006l ) and 
ikitzbichler fc White! 1 I2OOW) are indicated by a shaded region and 
a dotted line, respectively. 



(|2008[) . They used the GOODS-MUSIC catalog and their 
sample contains objects detected in z-band or i^s-band. 
The z-band and i^s-band limiting magnitudes are 26.0 
and 23.8 mag, respectively. The stellar mass function 
for our sample is in e xcellent agreement with that by 
lElsner et alT(|2008[ ). As lElsner et all (|2008l ) claimed, the 
stellar masses tend to be overestimated systematically if 
the IRAC data are not included in the SED fitting, espe- 
cially at z > 4. The discrepancy bet ween the stellar mass 
function of our sample and that by iDrorv et al.l ()2005f ) 
presumably attributes to this difference. Note that the 
redshift ranges of the resulting stell ar m ass functions 
of the sample by iDrory et al.l (|2005| ) and I Eisner et al.l 
(|2008( ) are from z = 4.0 to 5.0 and the representative red- 
shift (z = 4.5) is slightly lower than that {z = 4.8) of our 
sample. Also note that the discrepancy bet ween the stel- 
lar ma ss function of our sample and that bv lDrorv et al.l 
(|2005f) might be due to cosmic variance. 

The predictions of the theoretical models are also pre- 
sented in Figure fTHl We compare our fiducial stellar mass 
function with the pred i cted s tellar mass functions of the 
model by iBower et al.l (|2006D base d on the GALFORM 



Cole et al.l l 2000l) and the model bvlKitzbichler fc White! 
2007f ) based on ICroton et al.l (|2006D . both of which are 



semi-analytic models implemented on the Millennium 
Simulation and include the feedback effect from AGNs in 
the galaxy evolution. In the mass bin of ^ 10^^ Mq, 
our stellar mass function agrees with the theoretical mod- 
els. In the most massive part (M« ~ 10^^-^Mq), the 
number density of our sample is larger than the theo- 
retical predictions. However, if the mass functions of 
the models are convolved with Gaussian function with a 
standard deviation of 0.3 dex, which is a typical error in 
the SED fitting, by considering measurement errors, the 
models are matched with our result. In the most of the 
mass range (log(M*/M0) < 11.0), the number density of 
our sample is significantly lower than the models. 

5.4. The Stellar Mass Density at z ^ 5 

By integrating the derived stellar mass function, we 
can calculate the stellar mass density at z ~ 5. By inte- 
grating down to log(M*/M0) = 8.0, we obtained the stel- 
lar mass density of 7.0x IO^MqMpc"*^. As we discussed 
in §5.31 the choice of the assumed star formation history 
changes the shape of the stellar mass function, and hence 
the stellar mass density. In Figure [T71 we plot the stellar 
mass density calculated by integrating the fiducial stellar 
mass function (derived from the Category 1, 2, 3 objects 
with constant star formation history and excluding the 
overage objects). Assuming other star formation histo- 
ries makes the mass density larger: For instance, if the 
instantaneous burst model is assumed, the stellar mass 
density is 1.4 x IO^MqMpc"'^. Also, including the over- 
age objects makes the density larger by a factor of 2.7 
and 1.6 for the constant star formation model and the ex- 
ponentially declining model, respectively. However, the 
effect of including the Category 4 objects is not consid- 
erable: It makes the mass density larger by ~ 20%. If we 
restrict the maximum age of the models to 1.2 Gyr, the 
integrated stellar mass density is 1.1 x lO^MoMpc"^, 
which lies between our fiducial value and the plausible 
upper limit. 

The stellar mass density derived above is possibly still 
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Fig. 17. — Stellar mass density as a function of redshift. The stellar mass density from our sample is indicated by a filled star and 
the plausible upper limit is presented with a horizontal bar. Theoretical predictions from semi-analytical models are indicated by solid 
and dashed lines. The stellar mass densities of our sample, other observations at z < 4.5, and the theoretical models are calculated by 
integrating the stellar mass functions down to 10*Mq. 



underestimated because we miss LBGs with z' > 26.5 
mag and galaxies which are not selected by LBG se- 
lection. In ij5.31 we roughly estimated the contribu- 
tion from LBGs with z' > 26.5 mag to the stellar mass 
function. According to the rough correlation between 
UV magnitudes and stellar masses, the stellar masses of 
LBGs with z' > 26.5 mag are mostly less than~ 10^ AIq. 
By integrating the stellar m ass function estimated from 
the UVLF by llwata et all (2007) between IO^Mq and 
IO^ '^Mq, the contribution of the LBGs to the mass den- 
sity is calculated to be 4.6 x lO^MgMpc"^. 

We take the stellar mass density of 2.4 x lO'^AfoMpc^'^ 
derived with constant star formation models and includ- 
ing the overage objects and the contribution from the 
z' > 26.5 LBGs as a plausible upper limit. We show the 
upper limit indicated by a horizontal bar in Figure 1171 
The true stellar mass density at z 5 probably lies be- 
tween the value for our sample and the bar in the figure 

9 

We compare the derived stellar mass density from our 
sample to other observational studies . The stellar mass 
densities derived from the samples bv lGwvn fc Hartwickl 

® We estimated the uncertainty arising from the uncertainty of 
the stellar mass in the same way described in i|5.3l The average 
error in the stellar mass density is ~ 0.1 dex, which is negligible as 
compared with the uncertainty from the choice of star formation 
history. 



(120051). iFontana et all (I2006D . Pozzetti et all (120071) . 
iPerez-Gonzalez et al.l ()2008f). and Elsner et all (|2008l ) are 
obtained by integrating their stellar mass functions from 
lO^M© to 10"M pi. We also plot t he stellar mass den - 
sities at z ^ 5 bv lYan et al.l (l2006f). IStark et all (|2007f) . 
lEvles et all (|2007D . and lVerma et all (|2007D . These vaL 
ues are obtained by summing up all stellar masses (rang- 
ing from ~ 10 ^M(T) to ^ lO^^Mp ,) in their observations. 
The result by iGole et al] (|200lD is shown as the local 
value by integrating their stellar mass function down to 
10^ Mq. The estimation of the stellar mass depends heav- 
ily on the choice of IMF; the differing IMFs varies the 
mass systematically. In comparisons with other results, 
we applied corrections for consistency with our results 
in which we assumed the Salpeter IMF with lower and 
upper mass cutoffs of 0.1 and IOOM0, respectively. As 
illustrated in Figure I17[ our data point, including the 
plausible upper limit, is on a trend of gradual increase of 
stellar mass density with time. Our result agrees with 
other observations for galaxies at z ~ 5 jStark et al.l 
I2007t iVerma et all 120071: lElsner et all 120081 ) within the 
uncertainties. It is worth noting that the mass density 
at z ~ 5 is dominated by the massive end of stellar mass 
function. 

We also compare these observational results with 
th eoretical predictioii s of two semi-analytical models 
bv lBower et all (|2006l) and iKitzbichler fc Whitd (|2007l) . 
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The mass densities of these models are obtained by inte- 
grating their stellar mass function down to 10® M©. We 
applied corrections for the IMF as in the case of the com- 
parisons with other observations. Figure [17] shows that 
the number density for our sample is smaller by a factor 
of 2 — 3 than the semi-analytical models. However, if 
we take into account the plausible upper limit, our re- 
sult is almost comparable to the models. The theoretical 
models reproduce the overall trend of increase of mass 
density with time, though the number densities of the 
models tend to be somewhat larger than those of the 
observations in most of the redshift bins. 

The stellar mass density at z ~ 5 derived with a large 
sample of LBGs roughly agrees with that by the model 
predictions based on the CDM hierarchical structure for- 
mation scenario. However, since our fiducial stellar mass 
function disagrees with the model predictions as we dis- 
cussed in ^5.31 A larger amount of feedback to quench 
star formation might be needed in the lower-mass part, 
i.e., mass-dependent feedback process may be needed. 

6. CONCLUSIONS AND SUMMARY 

In this paper, we studied the stellar populations of 
Lyman Br eak Galaxie s at z ~ 5. We used the LEG 
sample by llwata et al.l (|2007f) obtained by the Suprime- 
Cam attached to the Subaru Telescope in the area of 
~ 500 arcmin^ around the GOODS-N, which consists 
of ~ 600 objects. For the mid-Infrared photometry, we 
used the publicly available data of the IRAC onboard the 
Spitzer in the GOODS-N. In addition, we observed the 
GOODS-N flanking fields (GOODS-FF) with the IRAC 
to cover the bulk (~ 80%) of the Subaru area. We se- 
lected ~ 100 objects which are isolated and not seriously 
contaminated by neighboring objects in the IRAC im- 
ages by eye inspection. For these objects, the rest-frame 
UV to o ptical SEDs were con structed. We used SEDfit 
package (jSawicki &: Yeel[l998t Sawicki, in prep.) to de- 
rive the properties of these galaxies: stellar mass, star 
formation age, color excess, and SFR. 

We assumed the constant star formation history, 
the metallicity of 0.2^0, the Salpeter IMF ranging 
from 0.1 M(7) t o lOO M^, and the extinction law by 
ICalzetti et al.l ()2000[ ), and found that the median val- 
ues of the stellar mass, age, color excess, and SFR are 
4.1 X lO^Mg, 25 Myr, 0.22 mag, and 141 MqYt-'^, re- 
spectively. The comparisons of the distributions of these 
pa rameters with those for t he z = 2 — 3 LEG sample 
bv lShaplev erall (|200lh and iShanlev eTall ([2001 . all of 
which are in the similar rest-frame UV and optical lu- 
minosity range, show the increase of the median stellar 
mass from z~5toz = 2 — 3bya factor of 3 — 4. The 
z ^ 5 LEGs are relatively younger by a factor of 10 — 20 
than the z = 2 — 3 LEGs. The median color excess of 
our sample might be slightly larger that that at z ~ 3 
and similar to that at z ^ 2. The median SFR of our 
sample is higher by a factor of 2 — 3 than in the z = 2 — 3 
LEGs. We suggest that the LEGs at z '--^ 5 are undergo- 
ing intense star formation making them dusty and they 
are dominated by younger stellar popul a tions than in the 
case of z = 2 — 3 LEGs. IVerma et al.l ([2007[ ) presented 
similar results with a smaller sample. If each LEG at 
z ~ 5 keeps the SFR derived from the SED fitting until 
z = 2 — 3, the expected distribution of stellar mass shifts 



toward larger than those derived at z = 2 — 3. This could 
imply that the SFR decreases from z^5toz = 2 — 3. 
We also compared the results for our sample with other 
studies for the z = 5 — 6 galaxies. Although we found 
similarities and differences in the distributions of the pa- 
rameters, we can not conclude their significance due to 
the small sample sizes of other studies. 

The large number of our sample galaxies allows us 
to derive the stellar mass function of LEGs at z 5 
after applying corrections for both z'-band and IRAC- 
band incompleteness. We compared the resulting stellar 
mass function with other observational studies. The stel- 
lar mass_Jmictionofour sample agrees with the result 
by iDrorv et al.l (|2005( ) in the most massive end. How- 
ever, in most of the mass range, the number densities 
of our samp l e are smaller than those of the sample by 
iDrorv et al.l ([20051). Meanwhile, our result agrees well 
with the result bv [Eisner et al.l ([2008[ ). The discrepancy 
bet ween the stel l ar ma ss function of our sample and that 
by [Drorv et al.l ([2005f l is considered to be due to the 
use of IRAC data in our analysis. We also compared 
our result with the predictions of semi-analytic models 
involving AGN feedback and found that although the 
number densities of our sample are comparable to the 
model predictions in the massive end of the stellar mass 
function, the observed number densities are smaller than 
those by the theoretical predictions in the lower-mass 
part. By integrating the stellar mass function down to 
lO*Af0, the stellar mass density at z 5 is calculated 
to be (0.7 - 2.4) x IO'^MqMpc"^. The stellar mass den- 
sity is dominated by the massive part of the stellar mass 
function. The stellar mass density of our sample is con- 
sistent with general trend of the increase of the stellar 
mass density with time obtained in other observational 
studies. Our stellar mass density is almost comparable 
to the models if we take into account the plausible up- 
per limit. The stellar mass density at z 5 derived 
with a large sample of LEGs roughly agrees with that 
by the model predictions based on the CDM hierarchical 
structure formation scenario. However, since our fiducial 
stellar mass function disagrees with the model predic- 
tions, some alterations may be needed for the theoretical 
models at high redshift. 
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APPENDIX 

A. EFFECTS OF VARYING STAR FORMATION HISTORY, METALLICITY, AND DUST 

EXTINCTION LAW 

In this paper, we adopted the models with constant star formation history, Z — 0.2Zq, and extinction law by 
iCalzetti et all (|2000t) . Effects of var ying the star formation history, metallici ty and dust extinction law have been 
explored for LBGs at lower redshifts (ISawicki fc Yeel[l998l : iPapovich et al.ll200lD and these effects, especially on stellar 
mass, are usually not large. Here we examine the effects of these various assumptions on the estimations of the stellar 
mass, age, color excess, and SFR for LBGs at z ~ 5. 

In the left sub-panel of Figure [T5k . the stellar masses derived by assuming the instantaneous burst and the expo- 
nentially declining star formation history with a time scale of r against the stellar mass derived for the case of the 
constant star formation history are plotted. To derive the stellar mass for the r-models, we used models with r=l 
Myr, 10 Myr, 100 Myr, 1 Gyr, and 10 Gyr. Although there seems to be no systematic difference (< 0.1 dex), the 
scatters are a — 0.59 and 0.35 dex for the instantaneous burst and the r-models, respectively. We also examine the 
two-component models, in which we put the additional star formation into the passively evolving component. We 
assume that the old component is the instantaneous burst model whose age is 890 Myr (z/ ~ 13). The spectrum of 
this old component multiplied by a flux ratio^° against the young component, is put into that of the constant star 
formation model. The stellar masses of some (~ 10) objects increase by ~ 1 dex, but those of most of the objects 
agree with the stellar masses in the case of constant star formation history with the median difference of 0.02 dex. 

We examined systematic effects of metallicity to the estimation of stellar mass. The middle sub-panel of Figure [T8k 
shows the output stellar masses with Z = I.OZq and 0.02^© against those with Z = 0.2Zq. Using solar metallicity 
models {Z = I.OZq) makes the stellar mass smaller systematically than using sub-solar models. Conversely, using the 
low metallicity models {Z = 0.02^©) derives systematically larger stellar masses. Mean offsets are —0.09 and 0.28 dex 
for Z = I.OZq and Z = 0.02Zq, respectively. 

We examined the effects of the choice of extinction laws on the stellar mass estimation. The extinction laws w e 
tested apart from the Calzetti law are SMC (jPrevot et al.lll984[ ). LMC, and MW extinction laws (|Fitzpatrick 1119861 ). 
These extinction laws are screen-type extinction, and we believe that the choice of the Calzetti extinction curve is 
appropriate, but we test these effects. In the right sub-panel of Figure [T5b. the derived stellar masses with using these 
three extinction laws are plotted against those with the Calzetti extinction law. The choice of these three extinction 
laws gives larger stellar mass especially in the massive part. On average, the SMC, LMC, and MW extinction laws 
give larger stellar mass by 0.07, 0.22, and 0.27 dex, respectively. 

Age, color excess, and star formation rate are also affected by the choice of star formation history, metallicity, and 
extinction law. We examined here these effects of varying the assumptions on the output parameters. 

As shown in the left sub-panel of Figure [T8b . varying star formation history affects the age estimation. On average, 
the age derived by assuming instantaneous burst, exponentially declining, and two-component models decreases by 
0.89, 0.70, and 0.44 dex, respectively, with respect to that by assuming constant star formation history. Most of the 
LBGs showing the old age (larger than 1.2 Gyr) in the SED fitting under constant star formation are younger than 1.2 
Gyr under other star formation histories. The middle sub-panel of Figure [T5b shows the effects of varying metallicity. 
Metallicity does not change the age estimation drastically as in the case of mass; on average, the age with l.OZ© and 
0.02Zo decreases by 0.22 dex and increases by 0.33 dex, respectively, as compared with the case of Z — 0.2Zq. The 
right sub-panel of Figure [T8b shows the effects of varying the choice of dust extinction law. The choice of other dust 
extinction laws generally increases the age systematically; on average, the age derived by assuming SMC, LMC, and 
MW extinction law increases by 0.48 dex, 0.40 dex, and 0.30 dex, respectively. 

Varying star formation history also affects the color excess estimation (left sub-panel of Figure fTSh). For most of the 
objects, the color excesses derived by assuming instantaneous burst are systematically larger than those by assuming 
constant star formation. There are two sequences; one shows no offset and the other shows a ~ 0.2 mag offset, and 
the average difference is ~ 0.1 mag. For some (~ 10) objects, the color excesses decrease toward E(B-V)~ mag. 
The color excesses by assuming exponentially declining star formation are also larger systematically by, on average, 
~ 0.05 mag than those by assuming constant star formation. There is no systematic difference (~ 0.03 dex) for the 
case of two-component models except for some (~ 10) outliers. The color excesses of some objects decrease for all 
star formation models, especially for the instantaneous burst models. The middle sub-panel of Figure [T8b shows the 
effects of varying metallicity. Metallicity does not change the color excess estimation drastically; there seems to be 
no systematic difference (< 0.01 mag) between the estimation of color excess by assuming I.OZq or O.O2Z0 and that 
by assuming 0.2^0. The right sub-panel of Figure [T8b shows the effects of varying the choice of dust extinction law. 
Assuming the SMC extinction law decreases the color excess systematically; though the difference is larger in the larger 
color excess, the average difference of the color excess from the Calzetti's law is 0.11 mag. On the contrary, assuming 
the MW extinction law increases the color excess; again though the difference is larger in the larger color excess, the 
average difference of the color excess from the Calzetti's law is 0.06 mag. Although assuming the LMC extinction law 
does not change color excess estimation systematically, there is a scatter of ±0.07 mag. 

Figure [TSH shows the effects on star formation rate. Varying star formation history (left sub-panel) also affects the 
star formation rate estimation but not so serious as compared with the age and the color excess, especially for the 

Note that r values of the exponentially declining models and the flux ratio in the two-component models are free parameters, and the 
degree of freedom decreases in the SED fitting. 
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Fig. 18. — Effects of varying star formation history, metallicity, and dust extinction law on the output parameters in the SED fitting. 
In panels (a), (b), (c), and (d), we show stellar mass, age, color excess, and star formation rate, respectively, derived by our fiducial model 
used in this work (abscissa) and by changing models (ordinate). In the left sub-panels, the output parameters obtained by assuming 
instantaneous burst (circles), exponentially declining (triangles), and two-component (squares) models are plotted against those obtained 
by assuming constant star formation models. In the middle sub-panels, the output parameters obtained by assuming I.OZq (circles) 
and 0.02^0 (triangles) metallicity are plotted against those derived with 0.2Zq metallicity models. In the right sub-panels, the output 



22 




Observed Wavelength (jim) Stellar Mass (witti Ha) 

Fig. 19. — (left): Best-fitted spectra of an objects with spectroscopic redshift with {solid line) and without {dotted line) Ha emission 
hne. The observed SED is indicated by filled circles. The SEDs of the best-fitted models are plotted by open triangles and open stars for 
models with and without Ha emission line, respectively, (right): Comparison of the best-fitted stellar mass with and without Ha emission 
line. 
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Fig. 20. — Effects of varying redshift on the output parameters in the SED fitting. In panels (a), (b), (c), and (d), we show stellar mass, 
age, color excess, and star formation rate, respectively, derived with fixed redshift (abscissa) and with free redshift (ordinate). 

larger star formation rate. The middle sub-panel of Figure [TSH shows the effects of varying metallicity. The effects 
on star formation rate is not large; the average difference between the star formation rate by assuming Z = I.OZq 
(0.02Zq) and that by assuming Z = 0.2Zq is 0.12 (—0.03) dex. The right sub-panel of Figure fTSH shows the effects 
of varying the choice of dust extinction law. Assuming the SMC and LMC extinction law decrease the star formation 
rate systematically (0.39 dex and 0.16 dex, respectively, on average). Although assuming the MW extinction law does 
not change star formation rate estimation systematically, there is a scatter of 0.33 dex on average. 

B. EFFECTS OF Ha EMISSION LINE 

About 70 percent of the SEDs of our sample LBGs show excesses in 3.6//m-band. This is very likely due to the 
presence of the Ha emission that shifts into the 3.6/im band. Thus we took into account the effect of Ha emission line 
in the SED fitting procedure. 

The Ha luminosity L(Ha) is calculated from the SFR of the given model by using the iKennicutt I (|1998f l relation: 

L(Ha) (ergs s"^) = 1.26 x 10''^ SFR (M© yr'^) (Bl) 

We examined the metallicity dependence of this relation. By using the BC03 models with metallicities of Z — 0.02, 0.2, 
and I.OZq, we derive the relation assuming all the ionizing photons are used for ionization and case B recombination. 
The systematic differences from the case of I.OZq are factors of 1.4 and 1.6 in the case of 0.2Zq and 0.02^0, respectively. 
For 0.2^0, the deviation from the Kennicutt relation is a factor of 1.3 and our estimation of L(Ha) is considered to 
be underestimated with this factor if the metallicity is 0.2Zq. Although the relation also depends on the assumed star 
formation history, the difference is less than a factor of 1.3 among r models. 

T he Ha lumino s ity ca lculated from its SFR is put into the model spectrum and the spectrum is attenuated with 
the ICalzetti et al.l (|2000[ ) extinction law. ICalzetti et al.l ()2000( ) argued the difference between the color excess for the 
stellar continuum Es{B — V) and that for the gas emission E{B — V), and presented the relation of Es{B — V) = 
{Q AA±Q .{)?,) E{B — V) . Thus, there is a possibility that the L(Ha) is overestimated by a factor of 2.3 in our prescription. 
While the L(Ha) may be underestimated by a factor of 1.3 for the sub-solar abundance, it may be overestimated by 
a factor of 2.3 due to the differing extinctions for the stellar continuum and the emission line. Considering these 
uncertainties, we adopt the Kennicutt relation. 
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Fig. 21. — Comparisons of the output parameters in the SED fitting for the z ^ 5 LBG sample bv lStark et al.l l|2007l ) with and without 
including NIR data (J and Ks). The photometric sample is indicated by circles and the spectroscopic sample is indicated by triangles. 
The comparisons of the stellar mass, age, color excess, and star formation rate are plotted in panel (a), (b), (c), and (d), respectively. The 
error bars are 90% confidence levels. 



In the left panel of Figure [191 show the best fitted spectra with and without H a emission line fo r an object 
(#144200). The redshift of this object was confirmed to be z = 4.69 by spectroscopy (|Ando et al.ll200j ). It is clear 
that the model with Ha emission fits the observed SED much better than that without Ha. We emphasize that we 
can obtain much better fit without increasing a free parameter. 

In the right panel of Figure [191 the comparison of the best-fitted stellar mass with and without Ha emission is 
presented. By taking into account Ha emission into the models, the stellar masses tend to be smaller than those 
estimated using models without Ha; there are two sequences, and the average difference is 0.32 dex. We examine the 
cases if Ha line is not included; the stellar mass function derived in i15.3l shifts toward larger mass systematically by 
~ 0.3 dex, and the stellar mass density discussed in §5.41 increases by a factor of ~ 2. If Ha emission is not included, 
the medina age decreases by 0.83 dex, the median color excess increases by 0.14 mag, and the median SFR increases 
by 1.00 dex, on average. 

C. EFFECTS OF VARYING REDSHIFT 

As we described in Sj3l we fixed the redshift of our sample objects to z = 4.8. Here we examine the effects on the 
stellar mass as well as other output parameters when we take redshift as a free parameter. We refit the observed SEDs 
of the Category 1 — 3 objects with the model SEDs, which are now parameterized by the redshift ranging from z = 0.0 
to 6.0 in a step of 0.1. Although we did not use V^-band photometry data in the fitting in the text, we incorporate 
y-band data here. Most of the objects (~ 89%) have the best-fitted redshifts within z = 4.6 ± 1.0, although some 
objects 9%) can be fitted as low-redshift {z ^ 2) objects like ellipticals as well. Eighty-four percent of the objects 
show no secondary minimum of value in lower redshift. These support that contaminants in our sample is small. 
If we exclude these possible low-redshift objects, the number density of the stellar mass function decreases by ~ 30% 
on average and the integrated stellar mass density decreases by ~ 30%. 

As illustrated in Figure I^Dk . there is no significance difference between the stellar masses with the fixed redshifts 
and those with the free redshifts except for some outliers. Excluding these outliers, the median of log(Af, (redshift 
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free))— log(M*(redshift fixed)) is —0.06 dex and the scatter is cr = 0.22 dex. Thus, even if we treat the redshift as a 
free parameter in the fitting, the change in estimated stellar masses is small. 

In the Figure [^Ob. I^Ob. andl^DH. the age, color excess, and star formation rate, respectively, derived by SED fitting 
taking redshift as a free parameter are plotted versus those obtained with fixed redshift of z = 4.8. Except for the 
outliers, fixing the redshift to z = 4.8 does not introduce a significant systematic offset and scatters are tr ~ 0.4 dex 
for the age and star formation rate and cr ~ 0.1 mag for the color excess. 

We additionally assess the uncertainty from the re dshift in the follow ing way: For each object, we assign the redshift 
randomly along the expected redshift distribution bv llwata et al.l (|2007t ) and run a large set of SED fitting. On average, 
the scatters of the obtained distributions of the parameters are 3 — 5 times smaller than the errors in the SED fitting. 

D. EFFECTS OF THE PRESENCE OF THE NEAR-INFRARED DATA 

In 21 we showed the results of the SED fitting analysis; the stellar mass, age, color excess, and star formation rate 
of 2 ~ 5 LBGs were derived. In this work, since we have no NIR data in the GOODS-N and the GOODS-FF, we used 
observed SEDs without NIR data. Thus we examined how large the discrepancy in the bes t fit parameters wi th and 
without N IR data. W e did t his test with the SED fit by using a sample of z - 5 LBGs bv IStark et all (|2007f l. The 
sample bv IStark et al.l (|2007f ) consists of 9 objects with spec-z and 34 objects with phot-z, for which both J and Ks 
data are available. 

Figure [m shows comparisons of the stellar mass, age, color excess, and star formation rate obtained with the J and 
Kg data and those without the J and Kg data. There seems to be no large difference between the stellar mass derived 
with and without NIR data; the stellar masses agree with each other within a factor of ~ 3. Thus, the estimation of 
the stellar mass is robust regardless of the presence of the NIR data. 

Figure [2Tb and [2lfc show the comparisons of the derived age and color excess, respectively, with and without the 
NIR data. Although the error bars are large, there seems to be no significant systematic difference between these 
derived parameters except for some outliers. We also plot the comparison of the star formation rate with and without 
the NIR data in Figure I^TH: the star formation rate seems to be rather securely determined except for some objects. 

However, the errors in the J and Kg data we use in the test are generally larger than those in other bands and 
the weights of the NIR data to the SED fitting are relatively small. The results of the test might be caused by this 
effects. Better constraints on the properties such as age and color excess are expected if sufficiently deep NIR data are 
available. 
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TABLE 1 

Photometry of the LBG sample (Category 1, 2, and 3) 
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.19 


± 


0. 


.11 


24.15 ±0.05 


24.37 ± 0.08 


133694 


27.56 ±0.12 


25. 


,21 ± 


0. 


.08 


25, 


.10 


± 


0. 


.07 


24.07 ± 0.05 


24.64 ±0.10 


134614 


< 28.09 


25. 


,69 ± 


0. 


.10 


25, 


.42 


± 


0. 


.09 


24.61 ±0.10 


23.99 ± 0.08 


135678 


< 28.09 


25. 


,20 ± 


0. 


.11 


25, 


.42 


± 


0. 


.12 


26.85 ± 0.62 


26.16 ± 0.47 


138613 


< 28.09 


26. 


,55 ± 


0. 


.21 


26, 


.36 


± 


0. 


.18 


< 25.93 


26.22 ± 0.49 


138763 


< 28.09 


25, 


,64 ± 


0. 


,10 


25, 


.69 


± 


0. 


,09 


23.81 ± 0.05 


24.72 ±0.14 


138810 


27.66 ±0.13 


26, 


,04 ± 


0, 


12 


25 


.88 


± 


0, 


11 


25.71 ±0.26 


25.44 ±0.27 


139906 


< 28.09 


26, 


,15 ± 


0, 


,20 


26, 


.09 


± 


0, 


16 


25.93 ±0.23 


9.99 ±9.61 


141088 


< 28.09 


25, 


,58 ± 


0, 


,09 


25, 


.22 


± 


0, 


,07 


24.62 ±0.10 


25.12 ±0.20 


141117 


< 28.09 


26, 


,52 ± 


0, 


17 


26, 


.31 


± 


0, 


12 


25.19 ±0.12 


25.53 ±0.21 


141368 


< 28.09 


25, 


,53 ± 


0, 


14 


25, 


.32 


± 


0, 


13 


24.82 ±0.12 


24.71 ±0.14 


142195 


27.89 ±0.16 


25, 


,20 ± 


0, 


10 


24, 


.95 


± 


0, 


10 


23.44 ± 0.04 


23.42 ±0.05 


144200 


27.02 ±0.07 


24, 


,15 ± 


0, 


,03 


24, 


.08 


± 


0, 


,03 


23.21 ±0.03 


23.68 ±0.06 


145330 


< 28.09 


26, 


,24 ± 


0, 


,15 


26, 


.13 


± 


0, 


14 


25.20 ±0.17 


25.18 ±0.21 


147153 


< 28.09 


26, 


.45 ± 


0, 


,23 


26, 


.27 


± 


0, 


18 


24.06 ±0.18 


< 24.14 


147965 


27.51 ±0.11 


25, 


.95 ± 


0, 


,12 


26, 


.08 


± 


0, 


13 


< 25.93 


25.76 ±0.34 


148198 


< 28.09 


24, 


.93 ± 


0, 


,05 


24, 


.50 


± 


0, 


,05 


22.47 ± 0.04 


22.66 ± 0.09 


149470 


28.16 ±0.21 


26, 


.50 ± 


0, 


,17 


26, 


.37 


± 


0, 


,12 


29.06 ±2.10 


< 25.64 


149604 


< 28.09 


26, 


,41 ± 


0, 


,19 


26, 


.43 


± 


0, 


,19 


< 25.93 


26.09 ± 0.44 


149667 


< 28.09 


26. 


,36 ± 


0. 


.17 


26, 


.36 


± 


0. 


.17 


23.98 ± 0.06 


23.48 ± 0.05 


152993 


< 28.09 


25. 


,87 ± 


0. 


.18 


25, 


.81 


± 


0. 


.10 


26.45 ± 0.35 


26.14 ±0.34 


156057 


27.61 ±0.12 


25, 


,31 ± 


0. 


.10 


25, 


.18 


± 


0. 


.07 


24.01 ± 0.06 


24.04 ± 0.08 


161503 


27.31 ± 0.09 


25. 


,70 ± 


0. 


.13 


26, 


.30 


± 


0. 


,12 


24.99 ±0.14 


< 25.64 


164830 


< 28.09 


25, 


,93 ± 


0. 


,15 


26, 


.07 


± 


0. 


,15 


24.41 ± 0.09 


24.77 ±0.15 



Upper limits are 3<t values in 1.6" diameter aperture. 

Errors arc Icr values. 

Upper limits are 3o" values in 2.4" diameter aperture. 



TABLE 2 
The best-fitted parameters 



ID 


logfM, (Mr,)]'' 


loff[Aee Tyr^l^ 

^^&L VJ /J 


E(B -V) (ma.s)'^ 


loelSFR fMovr~i')P 


A 


038819 




7 qq+2.39 
' •'^■^-1.14 


24+° °3 


2 40+0.40 


0.24 


038859 


^•-^^-0.79 


7 qq+l-f'f' 
' •'^■^-1.04 


"•^^-0.11 


i-75i;;;S 


2.20 


039340 


10-591!];;;;; 


9.78i;i;?;; 


n6+°'^-'^ 


0.951^;^^ 


0.01 


040064 




7-39l?i^ 


n 14+0.10 


1 54+0.69 


0.01 


046788 


10-94l°:tl 


9-57lOig 


nq+°'2^ 


l-O'J-0.37 


0.02 


048421 


10 14+"-"^ 


7 qi +2.18 
' -^-^-0.94 


n 9C+0.06 

"^•^^-0.25 


2 2q+0.42 
^•^''-1.16 


0.35 


048806 




7 qi +1.35 
' -^-^-0.83 


Q 1(1+0.06 


9 1 9+0.42 
^•-^^-0.76 


2.69 


050272 


io.27+;!-?|-^ 


7 4q+2.18 
' -^^-0.52 


n Q7+0.04 

^••^'-0.25 


2.8ll?;g 


0.11 


051334 


io.53i;;:| 


Q r)Q + 0.83 

"•^^-0.62 


n iq+o.o6 


9 00+0.34 
^■38-0.51 


9.78 


053312 


io.5oi;;;«^ 


8.84li;25 


'^■1^-0. 14 


1 75+0.75 
1- '^-0.61 


2.21 


061662 


Q 75+1-06 


6.76l?;22 


n qfi+018 
"^••^"-0.36 


2 qq+2 '14 


0.02 


062238 


11 lci+009 

11.10_g 77 


10.301°;°° 


n2+°-24 


1 01+1"2 


0.49 


063161 


7 qq+3'03 


c 79+4.58 
' ^-0.10 


n no+°'^^ 


Q Qy-j-0.38 

^■^ ' -1.59 


0.27 


063736 


11 45+"-i2 


7.6oi;;;^; 


o.74i;;;;;e 


oq+0.16 
'-'■"^-0.17 


1.07 


068263 


io.i7l!;;«l 




o.isi;;;;^ 


l-"^-0.79 


0.01 


070807 


ii-36i;;;;^I 


10.30+°°° 


0-07l°;23 


1 99+1.00 
-^•^^-0.12 


5.87 


071773 


in 62+""^'* 


9.68l°;i^ 


n nn+°-2o 

U.UU_Q gg 


1-^^-0.02 


1.21 


071987 


ii.07i;!;;51 


9.57l°;I? 


1R+0-50 


1 fi4+**'26 


0.02 


072556 


ii.28t;;;?;^ 


o 1 9+0.94 
°-l^-0.42 


c;i+0.06 


90+0.26 
■^■^■^-0.58 


0.17 


074486 


11.07l'^-il 


I0.20l'^;i° 


o.i2lg;^^ 




1.23 


074882 


10.03li|^ 


7.49l?;«i 


0.40l°;°*7 


9 c:7+0.44 
^•^'-1.56 


0.02 


075205 


11.211°;^? 


8.43lJ;?| 


0.38l°;°? 


2.86i;;;^2 


2.08 


077666 


lO-OllJi? 


6.66+?-«^ 


n 44+016 


Qc; + 2.25 
0.00_2 43 


0.01 


077684 




7.801?;;;° 


0.24l°;°9 


1 q2+0.70 
1-^^-1.14 


0.01 


078553 


9.62i;;f,?, 


6.45+f™ 


0.381°;°^ 


1 7+1.55 
■^■-'-'-2.47 


0.05 


079161 


9.60li;;^| 


7.28i?;;;| 


n Qq+0.06 
U.OO_Q 33 


2.34l?;« 


0.07 


079524 




8.64+;jt;4 


0.32l°;« 


1.93l?;?« 


0.01 


080726 


10.20+;!-:;? 


7-6oi;;i 


0.2812;?? 


2 64+0-25 


1.75 


082124 


ii-03i;;:;;^ 


7.08l'^;l° 


n 6q+0-02 


Q7+0.08 
-0.07 


2.02 


082855 


10.79l'^;^« 


q qq+0.31 


00+°'^* 

U-UU_Q gg 


q5+4.46 
"•'^°-0.02 


0.30 


083563 


9.50i;;g 


7.49l?;^i 


26+° °'^ 


2 n4+°'^'^ 


0.11 


083925 




7 4q+2-81 


n Qn+O O'' 
>J-5>J_0.30 


2.051?;!^ 


0.05 


084637 


10 n7+"'™ 

-^^•'J' -0.58 


8.53i;-56 


n 17+0.11 


i-62i;;;?? 


5.80 


084850 


9.68+;'';?; 


7.49+2:29 


r. 99 + 0.05 

'^■^^-0.22 


9 9-1 +0.50 
^■^-^-1.10 


0.47 


087802 


10.651;;;?;; 


7-9ilS:^t 


n 46+0.05 


9 on+o.2() 


3.59 


088084 


9-9313,;;;! 






1 58+l'°'' 
1-Oo_0.98 


0.01 


089484 


9.32i;;« 


7 ,Q+2.70 


0.24l°g;°| 


1 qfi+°-52 
i.yD_j 21 


0.12 


090105 




7 '^q+2-08 
'•■^^-1.25 


0.l8l°;?« 


9 1 5+0.47 
^■l'^-0.97 


0.03 


091420 


1 1 07+0.30 
-0.75 


q 78+0.52 
^- '''-1.98 


11 +0-25 
'^■11-0.08 


l-"'^-0.30 


0.31 


092240 


10 14+"'2l 


7 90+0.42 
'•^*'-0.31 


r, 99+0.02 
"■^^-0.02 


efi+0.16 
^■*»-0.19 


3.87 


092242 


^•°^-1.48 


7 oq+2.91 

7.39l°;^? 


2q+° °'' 

"•^"-0.29 


9 1 5+0.51 


0.06 


093559 


q 84+°'^^ 


2q+° °2 


2 48+°-13 
^•^°-().28 


0.42 


096484 


9-66l°;?f, 


7.28i;;?J 


rj 99 + 0.02 

"■^^-0.19 


2 oq+0.33 
^■•J^-1.02 


2.01 


096510 


11-681;;;;;^ 


I0.30l'>;°° 


f) 90+0.30 
^'■^'^-0.04 


-1 54+1-28 


0.09 


098022 


11-071°:°^ 


10.301°;°° 


n n2+"'"'^ 


n qq+0.25 
"■'"'-'-0.04 


3.97 


100184 


8.80l«;S 


90+1.46 


10+°°'* 


1 50+0.41 


1.88 


100509 


in 27+"-''" 

-^"•^'-0.58 


8.43ll;i^ 


0.i3l°g;i° 


1 q2+0-61 

-^•^^-0.59 


0.61 


102671 


q qq+o.ie 


7-6oi;;;^? 


0.l8l°;°2 


2 40+0. 16 
^■^'-'-0.23 


0.42 


103109 


11 11+0.31 
-^-^•-^-^-0.88 




n6+°'2'' 


1 47+1-26 
l-^'-0.25 


1.39 


103742 


9-43i;:?t 


7.39l2;»i 


97+0.06 
'^■^'-0.27 


2.071?;^! 


0.77 


103985 


7 q8+2'99 


c 79+4.58 
'^-0.10 


00+°'^^ 

U.UU_Q gg 


2-25l?;l 


0.28 


104189 


9-57i;;;?,j 


R q7+o.4/ 


-Ln no 
"••Ji_0.03 


2 fil +0-11 

^•"1-0.19 


0.57 


104766 


io.29i;;;?« 


7.80l°;J? 


30+° °* 


9 50+0.18 


2.38 


106944 


10 29+" '^'' 


8.64l?;«3 


oq+°'*'* 


1 74+0-79 

'^-0.40 


0.19 


107878 


in 81+0.10 
10.81_o og 


7.i8l'^:?J 


47+0.01 


(55+0.04 
■^■""-0.12 


1.03 


108167 


11 6S+°°* 

li.DCi_g gg 


10.301?;°° 


Q9-,+0.16 


1 49+°-65 
^•^"-0.08 


3.69 


108384 


11 14+°'^i 

^^•^^-0.23 


8 ni+i-*7 

8-Ul_o.52 


59+°-°* 


3.18+??! 


3.39 


108417 




00+4.47 
^•''•^-0.10 


o.oolS.f,^ 


2.29+?-'-^? 


0.21 


110593 


o c-n+0.71 

0.00_Q 29 


6 7fi+l-56 
D- 'D_0.42 


0.08l°;°^8 


-I 77+0:23 

'-0.78 


0.05 



29 



TABLE 2 — Continued 



ID 


log[M* (Mq)]'' 


log[Agc (yr)]=" 


i/(i3 — \/j (mag) 


log[SFR (M0yr ^jj^' 


2 

X 


112044 


10.79iOji 


9.571°;^^ 


o.i4l°o:?! 


1 oc + 1.58 
J.-'J'J—O.Sl 


0.01 


113060 


10.01+;;-?':? 


9.47l?-ii 


o-ool°:g3 


67+109 


1.37 


113120 


8.6ii»:g 


7 oq+l-35 


06+° °'' 

U.UD_g Qg 


1 9Pi+0.54 
J^-^'-'_0.47 


0.04 


113749 




7 r,Q + 0.00 

' -^O—O.OO 


0.001°:°° 


1 00+0.00 

1.00_Q QQ 


12.43 


115354 


10.4210-56 


1 0+1. 56 
"•^^—0.62 


0.39l°°l 


9 07+0.37 

•^••J' —0.94 


0.07 


115925 


»• J^'-'_0.27 


7.39l°;I^ 


(, f,r.+0.03 

U.UO_Q Qg 


1 7q+0.21 
J- '»_0.36 


5.79 


116678 




6.871°;^? 


r, 97+0.06 


•^•^^-0.34 


0.08 


116910 


9-971^52 




U.UD_Q Qg 


1 90+0. 78 
^^•^•3—0.33 


1.16 


117078 


9.83l»:i« 


7 90+0.31 

' •^°-0.31 


r, 09 + 0.03 

'-'•'J''_0.03 


9 C7+0.16 
''■J' —0.18 


0.24 


120190 


9.42lH'6 


6 q7+010 


iq+° °l 


2 46+° °^ 
^•''D_0.08 


2.35 


120554 




'•DU_2 50 


(, 24+0.25 
"-"•^^-0.24 


1.92+?-0« 


0.01 


120838 




7 70+0.52 


f, Qo+0.06 

U.UCi_g Q3 


1 26+°-3'' 

l-^D_Q 21 


15.24 


123533 


9.461^0^ 




n 29+°°'' 

>J-^c'_0.29 


11+1.44 
"^•^^—2.20 


0.79 


125665 


a-o' -0.12 


7.08l°;ll 


n 25+°°^ 


9 CI +0.13 
''•'JJ^_0.16 


0.73 


126010 


10.60l°i^ 


9.36l?;Ii 


(, (,7+0.18 

U-Uf _Q Q7 


1 oc+0.82 

1.00_Q 33 


0.27 


126510 


9.72lg:-2i 


^ qo+3.54 


0.42l°:°« 


7Q+0.96 
' ^-3.20 


0.02 


127245 


9.80tg:i? 


7 10+0.21 
'•-'-°-0.21 


0.28l°:°? 


2 64+°-" 
^■"^-0.08 


0.29 


127900 


8.86t°j« 


7 oQ+0.62 
' •'jy— 0.62 


n 02+°°^ 

'-'■U^_0.02 


1 c:n + 0.29 
1.0U_g 22 


7.08 


129670 


Q 70 + 0.39 

9- '0-0.88 


9.26+111 


'-'■'^'-'-0.00 


64+l l'5 

U.D^_Q ]^Q 


0.86 


130018 


Q 00+0.28 




i';+° °3 


2 16+°1* 
J^O_0.28 


2.40 


130851 


Q 90 + 0.62 


7 4q+l-56 
'•^y_0.52 


17+0-04 
'-'•l<_o.l7 


1 09+0.25 
J--O''_0.89 


0.94 


131482 


on+0.21 
22 


7 90+0.31 
' -^0-0.42 


n on+°°2 

U.UU_g QQ 


1 C7+0.26 
-0.08 


24.30 


133419 


9.68l«-J« 


7 70+1.56 
'•'U_o.42 


n iq+° °4 
>J-1»_0.19 


2 n'^+0-24 
z.uo_Q 92 


0.44 


133694 


9-36t°:25 


7.08l°1? 


n 2i+°°3 
"•^^-0.02 


2 29+°-l'^ 

^■^a_g j^g 


0.49 


134614 


10.46l°i^ 


9m+_tll 


00+°- 13 

U.UU_Q QQ 


r. Qo + 0.57 
'-'•^■3-0.04 


17.93 


135678 


8-49lH^ 


7 90+0.10 


00+° °° 

U.UU_Q QQ 


1 90+0.03 

±.^0_Q Q3 


7.02 


138613 


9-38t;:i 




ot:+0.09 
U.OJ_g 3g 


'^ n'^+1'5*^ 

O.UO_2 49 


0.05 


138763 


q 47+0.07 

y-^' -0.13 


r. -„ + 0.21 

6.76Io 10 


09 + 0.02 

'-'■>->''_0.04 


9 71 +0.06 
''^ ' J^_0.18 


1.14 


138810 


q 27+0.46 
a-^' -0.79 


°-^^-1.36 


n 07+°" 

U.U* _g Q7 


1 19+0.66 
^■^^—0.36 


2.38 


139906 


f;pr + 0.86 


7 oq+1.56 
'•■J»_1.04 


nq+0 06 
u.ua_Q Qg 


1 2q+0.51 
J^-^»_0.63 


0.17 


141088 


q 24+°-''2 
a-^^-0.29 


7 90+1.04 
'•^»_0.52 


r, 10+0. 03 


1.981^51 


7.19 


141117 


9.17±?,S 


7.39iJ;i 


"•^^-0.24 


1.801?:?^ 


0.42 


141368 


9.Sl+_lit 


8.74l?;g 


n 04+°-i* 


1 1C+0.73 

^•^°-0.23 


2.05 


142195 


10.38l°;« 


.39+1.35 

"•'->^_0.62 


(, 2S+0.07 
>J-^'J_0.18 


9 ,0+0.36 
''^ J^'->_0.79 


0.94 


144200 


9.8ll»:it 


7 90+0.21 
'•^°-0.21 


o.i9l°:°i 


Z.OO_Q j^3 


0.63 


145330 


q p;7+0.61 
-0.62 


Q QQ+1.35 


o.i4l°o:li 


1 01+0.66 

-^•''^-0.68 


0.22 


147153 




7.80l||° 


r, 07 + 0.08 

"■•->' —0.37 


9 90+0. 64 
^■^°-1.52 


0.01 


147965 


9.32l«;« 


8-641^;^^ 


n on+°'3'' 

U.UU_g QQ 


n 77+3.79 
'-0.03 


0.08 


148198 


10.471°;?^ 


7 oq+0.62 

'•■jy_o.3i 


n '!q+° °l 


in+O ll 
■^■^^'-'-0.33 


7.60 


149470 


8.031HS 


7 90+0.00 
'•^"-0.00 


00+° °° 

U.UU_Q QQ 


„ 77+0.00 

'-o.po 


9.72 


149604 


9.04i;it 


8-oii^;|? 


(, j(,+0.35 

>-'-lU_Q JQ 


^•^^-0.55 


0.01 


149667 


11.24tn^ 


10.301^5 


12+°" 
"•^^-0.03 


1.10±?,^«5 


4.01 


152993 


8.62l°:t? 


r- '7^1+0.52 
' • ''-'-0.83 


00+° °'^ 

U.UU_Q QQ 


o.96i°:?t 


2.24 


156057 


io.o7iS:g 


oQ~l~l-35 


1 7+0.08 
"J-^'-0.17 


1 09+0.46 
l-»^_0.77 


0.40 


161503 


np;+2.19 
o-UO_o.o4 


6.04t3;64 


o.oo1°q:J« 


2.02t°-20 


4.29 


164830 


Q ,-0+0.63 
».oo_Q 32 


7 

' ■'5y_0.62 


0.29l°:°? 


9 17+0:28 
^•^'-1.06 


1.11 



* Errors are 90% confidence and determined as follows: Monte Carlo realizations are computed for each object, and from the x'^ distribution 
of these realizations we determined the Ax^ which encloses 90% of the realizations. With this Ax^ , we define the 90% error contour for 
each object. Projections of the error contours onto the relevant parameter axes give the 90% errors on the individual parameters. 



